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1 Introduction 

1.1 Divergences of perturbation theory 

Very common situation which one encounters in physics is a lack of exact 
solutions. If Z{g) is some physical quantity considered as a function of a 
"coupling constant" [| g then the only what we have frequently at disposal are 
few first terms of perturbation theory, 

Z{g) = flo + aig + 025-^ + . . . + a„5'"- 

Unfortunately, without any regularization this expansion is badly defined 
in quantum field theory. In fact, in a formal perturbation expansion, the 
coefficients a„ are divergent and thus undefined. This is first kind of diver- 
gences which, in renormalizable quantum field theories^, can be cured by 
using a suitable renormalization prescription which gives a physically plau- 
sible recipe how to go from undefined unrenormalized coefficients an to the 
renormalized ones. We shall denote the renormalized coefficients by a^'s. 
This means that by renormalization one obtaines the perturbation series in 
a renormalized coupling constant g^, 

Z'^ign) = «f + afgn + a^gl + ... + a^gl, 

with well defined coefficients a^' s. In general, however, still one kind of diver- 
gences remains - it is the divergence of the perturbation series as whole. In- 
deed, if we do not consider the Yukawa model in two dimensions |jl| , fermionic 
field theories with UV cutoff where, due to the Pauli exclusion principle, 
the radius of convergence of perturbation theory for normalized (connected) 
Green's function is nonvanishing 0, as well as lattice gauge theories in the 
strong coupling region 0, perturbation theory of many other physical mod- 
els, especially without cutoffs and in the weak coupling region, is known to 
be divergent. The latter divergence is the main subject of the thesis. One has 
to know how to deal with such divergences, what is their origin, their phys- 
ical interpretation and significance if one wishes to have logically consistent 
understanding of quantum field theories. 

^Throughout the paper we shall call expansion parameter as a coupling. It may, how- 
ever, correspond to the inverse of physical coupling in the case of strong coupling expansion, 
to 1/A^ in the case of \/N expansion, etc. 

^including superrenormalizable ones 



In the case that a perturbation theory is asymptotic then it is no problem 
to deal with it in the region of small coupling (if we are only interested in 
perturbative phenomena). This is the case of asymptotically free theories at 
very large momenta or weak coupled theories at low momenta, respectively. 
If we are interested in nonperturbative phenomena then we stuck with the 
problem to give a meaning to divergent series. If the series is asymtotic (see 



section |1.3|) then by optimalizing the bound on the rest term (y) one gets the 
optimal number of terms to be kept. By keeping either more or less terms 
results in that our prediction based on the perturbation theory is getting 
worse. This optimalization of the bound on the rest term is for example 
essential to prove majority of our results. Indeed, our proofs follow mostly 
the same scheme and are in this sense a bit boring (sometimes were also 
boring for the author himself) : using the Euler-Maclaurin sum formula for 
a given sum, its asymptotic evaluation by the saddle point technique, and 
optimalization of the rest term. Nevertheless the scheme turns out to be 
quite efficient. 

In field theory the perturbation series cannot be truncated painless. When 
truncated perturbation theory becomes to depend on unphysical (renormal- 
ization scheme dependent) parameters. This freedom may be again used to 
optimalize the perturbative prediction and to extract a scheme independent 
result IQ, |[. Thus if a perturbation series is known to be asymptotic and 
divergent some time one can get a reliable prediction by keeping only finite 
number of its terms. However to establish the number of terms to be kept 
one needs to know the large order behaviour of perturbation theory. Infor- 
mation about large orders ceases to be only perturbative information and, 
in fact, it is a nonperturbative one. The fact that a presence of nonperturba- 
tive solutions can be seen by the large order behaviour of the perturbation 
theory could seem to be a bit surprising at a first moment but perturbation 
theory, if (whenever it is possible) suitably modified, is known to give even 
a good, nonperturbative approximation [^. As recently shown in |^ large 
orders of perturbation theory also determine multiparticle cross sections at 
asymptotically high energies. In order that this information be useful one 
needs a suitable summability method. This is the aim of this thesis to pro- 
vide a self consistent treatment of new result of the author on regular ana- 
lytic summability method which have been published in his several articles 

1,1,0, p. 

If one goes back to the history one finds that investigation of convergence 



of perturbation series in physical models was initiated by Dyson |]T2|, who 
considered the case of QED. More rigorous studies of the problem have been 
then done by a number of authors. The divergence of the pertubation series 
of the scalar bosonic \(j)\ model was proved by [|13|. More general polyno- 



mial bosonic field models and their correct large order behaviour have been 
studied, e.g., by |jl^, |T^. Perturbation theory for the Yukawa interaction 
(without cutoff's) in £ dimensions (2 < £ < 4) among fermions was shown to 
be divergent in [|16i . So far the divergence of the perturbation theory has been 



also proved for many quantities in other models of quantum mechanics and 
quantum field theory including the energy-levels of an arbitrary anharmonic 
oscillator [^^ |l^, the ground state energy of some nonpolynomial oscillator 
[|T^, one- loop effective lagrangian of QED [^, Green's functions in QCD 
[ pl[l and recently also some amplitudes of the bosonic string theory p2| . 

As for the explicite large order behaviour majority of results are based on 
the Lipatov argument [0], which gives in general the large order behaviour 
(for bosonic field theories) of the following type. 



an \< AaJ'rf 



'l + o{n)Yn\ 



(1) 



The bound was obtained by a formal steepest descent method in the Eu- 
clidean path integral formulation of quantum field theory for a variety of 
models |1^, |1^, |16|, |89|, 0. Unfortunately, the method is based on some 
additional assumptions which have not been proved [Q (cf. 1^)3) The re- 



sult, obtained by a phase-space expansion confirms the Lipatov analysis in 
< 4 dimensions, but in 4 dimensions renormalization was shown to disturb 
the Lipatov behaviour in such a way that the constant a = {3/2tt^) should 
be replaced by — /?i/2, where Pi{N) = {N + 8)/(27r^) is the first nonzero 
coefficient of the /3-function. Note that the ratio a/f3i tends to as A^ — > oo 
[ p5| , p^] . Rigorous large order behaviour of the massless (j)\ theory has been 



studied in [2/], but only for an UV-cutoffed theory. The rigorous large order 



behaviour has been only established in the case of a iV component massive 
A0| model. It is worthwhile to mention that the large order behaviour of 
bosonic string theory also differs from the Lipatov behaviour. Here the n-th 
order of perturbation theory is of magnitude {2n)\ in place of n\ in (bosonic) 
field theories. For further details see lEBI, |29|, |30|, plU and references therein. 



•^For quantum mechanical models as well as for models where horn-shaped singularities 
are absent works rather well. 



It is worth remembering that the factorial growth of perturbation theory was 
proved especially for bosonic theories while inclusion of fermions slow the di- 
vergence down (see, e.g., [0, |TB[)- Moreover, sometimes the factorial growth 



is an artefact of the approximation used and a more careful analysis may 
give even a convergent result [^. Similar divergences one also encounters 



in general relativity in connection with an existence of the de Witt integral 



Due to the field theoretical approach to critical phenomena |^^ initiated 
by Wilson |^ analogous divergences inevitably appear in the study of vari- 



ety of statistical systems. As for quantum spin systems with finite interaction 
the situation with regard to convergence is a bit better, at least in the high- 
temperature region which corresponds to the strong coupling region of lattice 
field theories [^ |36[. This is due to the fact that lattice models with finite 
interaction can be transformed on a polymer system |]3^ . Afterwords, by us- 
ing general results on polymer expansions I^Q], ^, one can derive convergent 
cluster expansions for them |P, ^, ^ . The convergence of cluster expansions 
allows to deduce upper and lower bounds on expectations of various types 
of observables such as Wilson loops, 't Hooft loops and others from bounds 
on polymer activities. The convergence also implies exponential decay of 
correlations as well as that above some temperature (coupling) there is no 
phase transition. This expansion is not, however, a power series expansion 
in general. In contrast to continuum field theory where one has general ar- 
guments for its large order behaviour |l^ there is no general argument how 



single terms of the cluster expansion should behave. In regard to convergent 



expansions in field theory and statistical systems see also [0, ^, ^ (and 
references therein). There a convergent perturbation expansion is obtained 
if, instead of the expansion in renormalized coupling constant, the expansion 
in powers of the running coupling constant is used. 

For complexity we note that some quantities of topological origin in su- 
perrenormalizable theories, like topological mass term in planar QED with 
the Chern-Simons term, can have finite perturbation series expansion with 
all coefficients but the first two being zero. This is ensured by the so-called 



non-renormalization theorem p4 |. 



1.2 Other perturbative techniquess 

In the last two decades a hard work has been done in developing perturbation 
theory where perturbation parameter is rank of the symmetry group of an 
underlying model. This approach was originally developed by 't Hooft [P5|]. 
He also noted a connection between the 1/A^ expansion and dual models 
[ ^ . If a connection were solidly established, many of the leading mysteries 
of QCD would be solved. The dual model has built into it confinement, with 
quarks at the end of a string. Also, a clear connection between QCD and 
the dual model would mean that the problem of dynamical mass generation 
had been solved, since the dual model certainly has a mass scale (the Regge 
slope). 

To elucidate the main features of the approach let us consider the familiar 
Hamiltonian of the hydrogen atom: 

2m r 

Since for e^ = we can solve this problem exactly - it is simply the problem 
of the motion of a free particle - one's first hope might be to try to understand 
the hydrogen atom in the small e^ regime by treating the interaction term, 
— e^/r, as a perturbation. This hope is frustrated because in the hydrogen 
atom e^ is not really a relevant parameter. It can be eliminated from the 
problem by redefining the scale of distances. 

After a rescaling r ^ tr, p —^ p/t, with t = 1/me'^, the Hamiltonian 
becomes 

H 

and one sees that the "coupling constant" e^ appears only in the overall 
factor me^ which serves merely to define the overall scale of energies and 
which could be absorbed in a rescaling of the time coordinate. Therefore, 
except for the overall scale of lengths and times, the physics of the hydrogen 
atom - and atomic and molecular physics in general - is independent of e^, 
and perturbation theory in e^ is meaningless. 

Atom and molecules can be described by the reduced Hamiltonian with 
e^ scaled out. For the hydrogen atom the reduced Hamiltonian is 



') 
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1" 
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The reduced Hamiltonian contains no evident free parameter. However, with- 
out a free parameter there is no perturbation expansion. Without a pertur- 
bation expansion what we can do? Suppose, since this is the case of QCD 
and of many other analogous problems that we were unable to diagonalize 
the Hamiltonian exactly, and then even a computer solution were imprac- 
tical. To make progress, we must make an expansion of some kind. Since 
there is no obvious expansion parameter we must find a hidden one. To find 
a hidden expansion parameter, we may treat as a free, variable parameter a 
quantity that one usually regards as given and fixed. 

For instance, we may take a cue from spectacular developments in the 
last years in critical phenomena. After decades in which the study of critical 
phenomena was frustrated by the absence of an expansion parameter, Wilson 
and Fisher suggested that to introduce an expansion parameter, one should 
regard the number of spatial dimensions not as a fixed number, three, but as 
a variable parameter |^. They showed that critical phenomena are simple 



in four dimensions and that in 4 — e dimensions critical phenomena can be 
understood by perturbation theory in e. Similar ideas were used by Bender 
at al. to study of A04 theory and some aspects of stochastic quantization 

How, by analogy, can we create an expansion parameter in our problem? 
Instead of studying atomic physics in three dimensions, where it possesses an 
0(3) rotation symmetry we may consider atomic physics in N dimensions, 
so that the symmetry is 0{N). One can show that atomic physics simplifies 
as N ^ oo and can be solved for large A^ by expansion in powers of 1/A^ 



4c 



Although this expansion in atomic physics is not very useful at A^ = 3 
(A^ must apparentely be at least six or seven for the 1/A^ expansion to give 
good result), illustrates the main features of the approach. In QCD, as in 
atomic physics, the coupling constant can be scaled out of the problem and 
the basic difficulty is the same as in atomic physics - the seeming absence of 
an expansion parameter. For further details see Q . 



So far the 1/A^ expansion has been successfully used to solve the Gross- 
Neveu model [0], in the study of two-dimensional Yang-Mills theory [^] , 
f/(oo) lattice gauge theory ^|, the generalized two-dimensional U{N) Thir- 



ring model [^, and the 0(A^)-symmetric sigma model |^. In last years re- 
sults of 1/A^ expansion have been also succesfully applied to two dimensional 



quantum gravity (see \E^ and references therein), as well as to quantum an- 



tiferroinagnets in connection with high-temperature superconductivity [^ . 



As for the study of planar field theories (A^ -^ oo limit) see |57, ^. The 
1/A^ expansion can be improved by a suitable summability method just as 
conventional perturbation theory [^ . 

Unfortunately, for Yang-Mills theory, the N = oo limit does not seem to 
be exactly soluble, even though it is significantly simpler. Nonperturbative 
approches to Yang-Mills theory are sorely needed in particle physics. There 
is not much hope that the Yang-Mills theory itself is exactly soluble, so one 
has to look for another theory which is "close enough" to it. Then Yang- 
Mills theory can be studied by an expansion around this soluble theory. As 
a general rule the more symmetry a theory has the more easily it can be 
solved. 

Can we enlarge the symmetry group futher and thus find a soluble theory? 
The answer is "yes" . Consider Yang-Mills theory with structure group G on 
a spacetime X x R. In the gauge Ao = the structure group is the set XG 
of all smooth maps from X to G. Let us consider G = U{N). When we 
take the N —>■ oo limit in the usual way, we consider an infinite sequence of 
Yang- Mills theories with structure groups 

U{1) C f/(2) C . . . U{N) C U{N + 1) C . . . U{oo). 

The gauge group of Yang-Mills theory in the N = oo limit is then XU{oo), 
which contains XU{N) for any A^. But there exists a group (the universal 
gauge group) f/2 which contains XU{N) not only for all A^, but also for all 

There is too early to say whether this approach will be successful, but 
it by no means provides a useful step to find an exactly soluble theory for 
perturbation theory calculations. 

1.3 Asymptotic expansions and strong asymptotic con- 
ditions 

As it was discussed above perturbation series in quantum theory are mostly 
divergent and can have at best the meaning of asymptotic series. Asymp- 
toticity means in general that the perturbation theory cannot determine the 
solution uniquely. More precisely given an arbitrary sequence {a„}J^o of 
complex numbers and an arbitrary sector- like domain D, there exists for 



some e > function f{z), which is regular in D,. := D n {z \\ z \< e} and 
Urn {f{z) -ao- ...- a„2;")/z"+^ = a„+i 

exists, or equivalently 

oo 

f{z)r^J2^nz'' iz^O,zeD,), (2) 



i-G-) J2'^=o(^nZ^ is an asymptotic series of f{z) in the region D^ [^ |6^. In 
general there are infinitely many functions with the above properties. 

Asymptotic series may be divergent as well as convergent (note that a con- 
vergent asymptotic expansion does not prevent a function from being singular 
at the origin) . To deal with such series one has to look for a maximal region 
D of the complex coupling constant plane in which the asymptotic expansion 
is uniform. In order that such series determine its sum f{z) uniquely this 
series has to satisfy some additional conditions, so called strong asymptotic 
conditions (SAC). SAC are then relations which put into the balance the 
shape of the region D, or analyticity in the complex coupling constant plane, 
with some bound on the rest term Rn{z), 



N-l 



f{z)= Y^a^z^ + R^iz), (3) 

n=o 

in such a way that only one function f{z) can satisfy them. SAC are con- 
ditions which by definition prevent appearence of nonperturbative terms like 
exp{—A/g) in perturbation theory (or a typical "tunnelling" like amplitudes), 
i.e., terms whose perturbation expansion is identically zero. Indeed, to be 
unseen in the standard perturbation theory the (non-zero) contribution of 
such a term should have an asymptotic series which is identically zero. Such 
are for instance all nontrivial minima A of action S for which S{A) < cxd. 



like instantons in non-abelian theories |Q or kinks in the double-well anhar- 



monic potential [^, |BB[, etc. Violation of SAC then indicates a presence of 
nonperturbative effects in a theory and instability of its ground state. It was 
emphasized that the violation of SAC is a more serious problem of a theory 
than the divergence itself, therefore the violation of SAC leads to ambiquity 
of perturbation theory predictions fll], |6^. Of course, this ambiguity de- 



pends on the behaviour of the rest term Rn{z) in (^ for N — »• cxd. As it has 



been exposed above, optiinalization of the bound may determine an optimal 
number No of terms to be kept and the value of Rn^{z) then determines a 
maximal error by which the exact solution differs from the perturbative one 
even in the case that the series violates SAC and does not determine a unique 
solution. Sometimes ambiguity of the perturbation theory predictions, e.g., 
in bosonic string model are highly desirable, because "if perturbative string 
theory were make sense, string theory would have nothing to do with phys- 



ical reality" |2^. This is because there are many features of perturbative 
treatment of string theory which are not shared by the real world. When 
SAC are fulfilled then the coefficients a„ of perturbation theory determine 
f{z) uniquely. To obtain this function an appropriate analytic regular sum- 
mability method can be used p2|. We say that a summability method a is 
regular if and only if 

oo oo 

a(5:a„z") = ^a„z", (4) 

n=o n=o 

whenever the r.h.s. of (^ converges. Analogously, we say that a summability 
method a is analytic if for every power series X) o-nZ^ with nonzero radius of 
convergence 

oo 

a(Ea„^") = /(z), (5) 

n=o 

whenever the l.h.s. of (^ exists {f{z) now being an analytic continuation of 
X]a„z"). The moment constant summability methods as defined by |B^ are 
automatically regular, although not analytic in general. We confine ourselves 
to the class of analytic moment constant summability method (AMCSM) 
|6^ , one member of which is the well-known Borel summability method, 
frequently used in physics (see for a recent review). 

To illustrate the main features of the AMCSM let us consider a very 
simple example of the series I].^o z"^. Using the identity 

/ e-'e dt = 1, 

n\ Jo 
it is always true that 

oo oo n „oo 



within the radius of convergence of I]^o-z", i.e., within the disc \\z\\ < 1. 
However within the disc of convergence one can always interchange summa- 
tion and integration. Thus, 

E-" = / e-\Y.^-^)dt. (6) 

n=o •^° n=o "'■ 

Here, the r.h.s. of (|) is nothing but the Borel sum of J2T=o^"'- I^ this case 
it can be calculated explicitely, 

e-*l^ dt= e-'e^' dt = e-'^'-'^ dt = ^, (7) 

n] Jo Jo 1 — z 



o 



whenever Rez < 1. Thus, we have obtained the expression which converges 
ioT Rez < 1, i.e., outside the disc of convergence of I]J^o^" ^^^ gives there 
an analytic continuation of Y.'^=o ^" (^^e Fig. ^. 

1.4 The Borel summability method 

1.4.1 The Nevanlinna theorem 

A main mathematical tool for investigation whether a series is Borel sum- 
mable or not provides the Nevanlinna theorem |^, |68| : 



Theorem 1 .■ Let f{z) be analytic in the circle Cr := {z \ Rel/z > l/R}, 
continuous up to the boundary, and satisfy there the estimates 

f{z) = Y, akZ^ + Rn{z), 

k=o 

with 

\RNiz)\ < Aa^ m \ z l"" , (8) 

uniformly in N and in z & Cr. Then B{t), 

oo j.n 



10 



converges for | t |< 1/cr and has an analytic continuation to the striplike 
region S^ '■= {t \ dist(t,R^) < 1/cr}, satisfying the bound 

\ B{t) \< Ke\'\/^ (9) 

uniformly in every S^' with a' > a. Furthermore, f can be represented by the 
absolutely convergent integral 

f{z) = {l/z) e-"^B{t)dt (10) 

Jo 

for any z G Cr. 

Conversely, if B{t) is a function analytic in S^" (cr" < cr) and satisfying 
there the bound ^, the function f{z) defined by j^Wj ) is analytic in Cr and 
satisfies (^ [with an = B^'^'it) \t=o] uniformly in every Cri with R' < R. 

1.4.2 Application to physical models 

The reason why the Borel transform has been so successfully used in physics 
is due to the fact that the path integral can be rewritten in the Borel form 
[^, ^. Indeed, consider a euclidean functional integral which has been 
rescaled as in a loop expansion 

Cig')=J[dA]eM-SiA)/g% 

Functional integrals are the N ^ oo limit of integrals like 

N 

g-^'XjliJ dxiexp[-S{x,)/g% 



where Af is the usual normalizing factor. Using the Cauchy integral formula 
this can be rewritten as 

dz ,12^ r(A^/2 + 1) 



AfjY[dxij> —exp(-z/g 



[S{X,) - 2]^/2+l ' 



where C is a contour encircling the range of S, which is normalized to run 
from to oo. An interchange of integrations is allowable for finite N and 



11 



reasonable actions yielding 

G{g') = / ^exp(-^//)i?;vW, 



^^ ,^ r(iv/2 + 1) 



B.i.) ^ Nj nf.„c/x.^^A_^-p^ . (11) 

A final step is to take functional integration limit, B{z) = liniAr_^oo-BAf(^)- 
Note that since denominator is raised to power N/2 + 1, naive power counting 
guarantees convergence even when S is only quadratic (as in free theory). 
On the other hand the coefficients a„ of the perturbation series in quantum 
models do not grow faster than A(T^r{an + /5), where F is the usual gamma 
function and the constants A, a, a and /? do not depend on n, what is just 
the case which can be dealt with the Borel method. 

The Borel summability of the perturbation series of the physical quan- 
tities has been studied by many authors. In quantum theory it was firstly 
used by ||70| to show that one-loop effective lagrangian of QED exactly cal- 
culated by Schwinger in the proper-time formalism [|7T| is for the vanishing 
electric field nothing but the Borel sum of the perturbation series (see also 
7^). The Borel summability of the perturbation series has been also proved 



in the case of the energy-levels of an anharmonic oscillator in any finite di- 
mension [jl^ [T^, |73|, the resonances of the Stark effect on a Hydrogen-like 



atom |]72|, 74 ] , in the case of the Zeemann effect [^ , for some nonpolynomial 



oscillators [0 and for the planar asymptotically free massive field theory 



for sufficiently small coupling ||5^. The method was also used for a rigorous 
perturbative construction of the massive euclidean Gross- Neveu model [Q]. 
In combination with conformal mapping or with the Fade approximation the 
Borel summability has been successfully used to compute the critical coeffi- 
cients of the statistical models |7^, the zeros of the beta function [^, the 



ground state energy of the anharmonic oscillator at infinite coupling limit 
7^ , to compute the beta function of the massive A0| model [^ or to sum 



dominant coupling constant logarithms in QED-^ [|y], etc. 

1.4.3 Structure of the Borel transform 

In physically realistic models the structure of the Borel transform is very 
involved. The obstacles for the Borel summing of a perturbation series ap- 
pear as singularities of the Borel transform on the real positive axis. These 

12 



properties of the Borel transform can be deduced as follows. The integral for 
Bn{z) is very much like a Feynman parameter integral. Thus the Landau 
conditions may be applied to ask when Bi\f{z) is singular [^. The Landau 



conditions have a direct interpretation in terms of solutions to the euclidean 
equations of motion: Bj\f{z) may be singular for z such that there exist x^ 
for which 

z = S{xt), 

= dS/dx, U^=,, . (12) 

This means B]sf{z) may be singular for z equal to the action of a solution 
to the classical euclidean equations of motion (discretized). This is true for 
complex as well as real-valued solutions. Such solutions are generally called 
pseudopartides or instantons. It is widely supposed that in the large- A^ limit 
the solutions to the discretized equations approach a limit which is identified 
with the solution to the continuum equations of motion. 

In general, attached to each singularity will be a branch cut. Because 
B{z) is integrated along a contour C encircling the positive real axis, it is 
convenient for real, positive singular points to have their branch cuts running 
along the positive real axis. 

The only singularities on the first sheet are on the positive real axis. 
These can arise from real solutions with real positive action. When no in- 
stantons are present the only singularity on the non-negative real z axis is 
at 2; = 0, the singularity arising from the vacuum solution. There is a cut 
extending from this singularity to positive infinity. When instantons occur, 
there will be another contributions to the Borel sum othe than those visible 
in perturbations about vacuum. The instanton singularities will have branch 
cuts attached to them. The total discontinuity of B{z) will be a sum of 
discontinuities across branch cuts. 

In renormalizable (not superrenormalizable) theory there are rather per- 
suasing arguments that apart from instantons another singularities can ap- 



pear which go under generic name of renormalon singularities [21, 82, S3|. A 



toy, finite dimensional model where renormalon singularities appear has been 



investigated in p4| in the context of Wilson's renormalization group acting 
in the space of hamiltonians of the cutoff system. The renormalon singular- 
ieties have been shown to be related to the presence of a marginal direction 
at the fixed point of the renormalization group besides the relevant or irrel- 
evant ones, or to the resonance in the dynamical system language. Despite 

13 



Figure 1: Supposed structure of the Borel transform of QED. Here the units 
are Stt if a is the original expansion parameter. 

great effort presence of renormalons in renormalizable quantum field theories 
is not rigorously established yet p5| , |26[| . This important unsolved problem 



is in a close connections to the so-called triviality problem - all known four- 
dimensional renormalizable quantum field theories are expected to be trivial 



pq , pmi 



The suggested structure of the Borel transform of QED and QCD is 
shown in Fig. |l| and ^ Note that the position of renormalon depends on 
the sign of the first nonvanishing coefficient of the corresponding /5-function. 
While in an asymptotic free theory UV-renormalons are harmless (lie on 
the negative real axis) in nonasymptotic free theory like QED renormalons 
prevent the Borel summability (lie on the positive real axis). As for QED 
see also |^5[. (We would like to stress out that the calculations which led 



to the renormalon [|2T], ^ was only performed for specific subclass of all 
graphs, which contribute to the given order of a theory, and may only serve 
as indication but not at all as a proof). 

1.4.4 Shortcomings of the Borel summability method 

Unfortunately, the use of the Borel summability method is sometimes very 
limited. Nonsummability means that perturbation theory violates strong 
asymptotic conditions (SAC), or, what is equivalent that perturbation the- 
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Figure 2: Borel z 'plane for QCD. The circles denote IR divergences that 
might vanish or become unimportant in colour-free channels. 

ory predictions can not be unambigious. This happens due to presence of 
nonperturbative effects which cause perturbation theory around the trivial 
minimum to be unrehable. Nonperturbative effects do appear and destabihze 
a trial (classical) ground state in models having classically degenerate ground 
state IQ which is the case, e.g., of the double well potential |^, ^ ^, gauge 



field theories |^4|, |6^, ^, as well as of the heterotic string and superstring 
[Q. However, SAC can be also violated in the case of perturbation the- 
ory around a stable nonperturbative ground state (after the nonperturbative 
effects have removed the degeneracy of the ground state [^), as has been 
shown in the case of the double- well potential |[87[|P|. The physical meaning of 
the non Borel summability is illustrated on some simple models in [^, |91 . 



The next shortcomings of the Borel summability method is connected 
with its analytic properties. Indeed whenever the Borel sum exists then it 
defines a function regular in some sector-like neighbourhood of the origin (the 
convergence of the Borel integral need not be absolute!) |9^. On the other 
hand a maximal region of analyticity of realistic field theories in the coupling 
constant plane may be very small. For instance, there are strong arguments 
that two-point Green's functions of four dimensional renormalizable mass- 
less field theories are only analytic in a horn-shaped region: a wedge bisected 

^I thank J. Magnen for acquinting me with this reference. 
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Figure 3: Expected form of the maximal region of analyticity in the coupling 
constant complex plane of four dimensional renormalizable field theories. 

by the real axis and bounded above and below by circles which are tangent 



to the real axis at the origin ||21], ^ |9J] (see Fig. y). Hence the maximal 
region of analyticity has a cusp, i.e., zero opening angle at the origin. The 
appearence of a horn-shaped region of analyticity is in fact a nonperturbative 
result which we have about a theory. It is not seen in a perturbation theory 
and arises if the analytic structure of the (response) Green functions for com- 
plex momenta p^ is combined with asymptotic freedom. This is the essence 
of 't Hooft's argument pT| , ^^ which relates the standard momentum space 



branch points for timelike p^ > to singularities in the complex plane of a 
suitably chosen coupling constant by solving the Callan-Symanzik equation 
(for generalizations to more general couplings see [^). In fact, the use of the 
Callan-Symanzik equation is not crucial in 't Hooft's argument and it goes 
through even if a theoretical evidence of the asymptotic freedom is replaced 
by an experimental evidence (if available) of the asymptotic freedom, i.e., 
that a measurable coupling g'^{fi) behaves like 

^9'{f^) = -(3i9\f^) + Oig%f^)), (13) 
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on a sufficiently high mass scale /i. Note that in this case we always can write 
an analogue of the equation ([T3|). Provided that such information is available 
one gets rid of the problem of definition of a suitable coupling and the horn- 
shaped region appears as UV effect once one assumes that at high momenta, 
measurable quantities are only functions of the ratio p^/^l, Ho being an 
integration constant of (|^). Analytic structure of response functions is by 
itself a nonperturbative information about a theory which arises when very 
general theoretical principles of unitarity and causality are combined with 
experimental data. Hence, the appearance of horn-shaped regions is in fact a 
very general phenomenon resulting whenever the following points are satisfied 
or available: 

1) unitarity 

2) causality 

3) asymptotic freedom (either theoretical or experimental evidence) 

4) experimental data (which determine the position of thresholds and branch 
points in the momentum plane). 

Some other arguments for producing singularities in the ^f-plane are not at 
all related to branch points in the momentum plane, and are due to the 
separatrix of the differential equation for the running coupling constant, i.e., 
to the line in the gf-plane which separates the asymptotically free region 
of initial values g from the non-asymptotically one [p6| , [97[] . If the Callan- 
Symanzik procedure for calculating the asymptotic behaviour of one-particle- 



irreducible Green's functions in the deep Euclidean region ||9^ is correct 
then such horn-shaped regions of analyticity is a typical UV effect of all 
four dimensional renormalizable field theories. What is interesting is that a 
horn-shaped region is also the maximal region of analyticity in the complex 
coupling constant plane of energy levels of a simple quantum mechanical 
model - the anharmonic oscillator in the massless limit |9^ . 



If one imposes UV cutoff then two- and four-point Green's functions of a 
four dimensional massless field theory can be analytic in a disc as was exem- 
plified in a rigorous construction of the critical (f)\ theory [^. However, if one 
removes this UV cutoff the disc should shrink to zero. Recently it was shown 
that some other quantities, like /5-function (defined perturbatively using the 
BPHZ renormalization with subtractions at vanishing external momentum) 
of the massive 04 model without UV cutoff may be Borel summable and 
hence analytic in a disc |7^. For the sake of completness we recall that in 
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the case of 0^ with d < 4 there is no horn of singularities [^ . As for the non- 
gauge field theories some arguments in favour of ambiguity of renormalized 



04 field theory has been given due to 1/A^ expansion [^ . 

The Lipatov bound (|ID is compatible with SAC only for a regions of 
the asymptotic type (0, rj) (a is the same for all quantities in a theory 
11^ , pTS] , p!6| , [89| , ^^). As we have mentioned above a general rigorous jus- 
tification of the Lipatov method has not been established so far. However 
further investigations would not change the fact that perturbation theory of 
four dimensional renormalizable field theories violates SAC. If the large or- 
der behaviour of the model in the limit of vanishing UV-cutoff will be proved 
to violate (^) then this will allow for nonperturbative solutions to the the- 
ory and can sheds new light on the triviality problem of the above theories 
|10CI| , |101|| . To get reliable perturbation theory one must find a nonperturba- 
tive ground state and expand around it ||1(J2|| . In the case of the double well 
potential also some other techniquess were advocated ||103|| . However, we 
cannot exclude the situation that the coefficients a^ of perturbation theory 
will obey the bound (pB]) in this model. Provided it is so than no nonper- 
turbative solutions can occur and the Borel transform should be an entire 
function, what contradicts expected singular structure of the Borel trans- 
form. However, it may well happen that there is no reason to look at the 
Borel transform of the four dimensional renormalizable massless field theory, 
because in this theory convergence of the Borel transform contradicts the 
observed horn-shaped region of analyticity [P, ^, in distinction from the 
lower dimensional case 1961 



2 Goals and methods of the thesis 

2.1 Goals 

The main goal of the thesis is to overcome some shortcomings of the Borel 
summabihty method. Let us summarize the shortcomings which have been 
solved by the author. 

A) the Borel method cannot cope with the horn-shaped singularity, exhibited 
in the massless QCD [pl|], in the "massless" limit of the anharmonic oscilla- 



tor |^4[, and in some other models. (For recent status of the problem see ^6 
and references therein.) 



Another shortcoming of the Borel method is the following. 

B) In the regular case, i.e., when the series on the r.h.s. of ([^ has a nonzero 
radius of convergence, it can in general be analytically continued onto a re- 
gion of the complex plane which is larger than the region of z for which the 
Borel sum exists. If the analytic continuation f{z) is not an entire function, 
the Borel sum may not exist for all z from the Mittag-Leffter (principal) star 
of f(z) [^] (further MLS{f)), because always some sector-like domain has 
to be discarded from the complex plane 0, |8ll . 



To illustrate point (B) note that the standard Borel method sums the 
series I]^o ^" only in the complex halfplane Re z < 1 but the analytic con- 
tinuation 1/(1 — z) of the series exists in the whole complex plane except 
the point z = 1. We should like to demonstrate in simple examples that the 
Borel summabihty method does not have any priviliged role and that it is 
just an ordinary method among a continuous variety of other analytic mo- 
ment constant summability methods (AMCSM). It should be used when a 
divergent series has to be summed up to a function regular in a disc K{0, R). 
Provided the sum of a given divergent series is known to be analytic in some 
other region then the corresponding AMCSM compatible with this analyticity 
requirement has to be used. 

To generalize known SAC to horn-shaped regions which may be of some 
physical importance is another goal of this thesis. The well-known example 
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Figure 4: Region of the Borel summahility of 1/(1 



of SAC provides the Nevanlinna theorem [^, ^. These conditions are, 
however, apphcable only when sum of a given perturbation series is to be 
analytic in a disc /^(O, R) := {z \ Re{l/z) > 1/R} tangent to the imaginary 
axis at the origin (modulo mapping z -^ z^'"'). On the other hand, we have 
seen that many theories are known to be analytic in a horn-shaped region 
only. Without the validity of SAC one cannot prove that a statement which 
is valid order by order in the perturbation theory is a property of the full 
theory. One such example is, e.g., the Callan-Symanzik procedure where 
certain mass insertion terms are ignored and assumed small even though one 
has been able to prove this fact order by order. It was found ||93| that the 
Borel summability or SAC given in |^, ^, ^ cannot help in justifying the 
Callan-Symanzik procedure just because of the horn-shaped singularities. In 
[ P^ it was stated that the horn-shaped singularity puts the Callan-Symanzik 
assumption about the mass insertion term in doubt. Our point of view is 
that there is no controversy between the horn-shaped singularity and the 
Callan-Symanzik procedure since one has to use SAC appropriate to a given 
analyticity region. 

We shall also give a further improvement of the SAC as discussed in 



20 



[|6^ , |67| , |68| , and show that a product of functions obeying these SAC also 
obeys the same SAC. The knowledge about the relation of SAC and sum- 
mability of perturbation theory is important in order to understand, e.g., 



the recent result on divergency of the bosonic string perturbation theory p2 



and frequently quoted in connection with two dimensional quantum grav- 
ity since there no care is taken about analyticity to which the perturbation 
theory should be summed up, and physical consequences are simply drawn 
from Borel nonsummability. Similar "negligence" of analytic properties also 
occurs, e.g., in [|l^, ^ or |2^. Fortunately, the series treated in ||15|, ^ ^| 



can be dealt by Lemmas ^ and |^ (see below). As for |2^ the results of 
this paper can be valid provided a fixed-angle elastic scattering amplitude 
of bosonic string A{s,(j);g) (see also [^]) is analytic in a region D of the 
sheeted complex plane with an opening angle B > Qtt since the resummed 
amplitude Aresum is such [|rB], |57[ . 



SAC discussed below may provide a tool to decide whether such ambiguity 
(degeneracy of the ground state of an underlying theory) is actually present. 
However, some futher work has to be accomplished in order to gain the 
maximal domain of the analyticity, establish a rigorous large order behaviour 
or to give a suitable lower bound on the coefficients of the perturbation 
theory. 

2.2 Methods 

To solve the above shortcomings (A) and (B) of the Borel method we use 



theory of complex variables [^, asymptotic series ||104| , |105|| , and divergent 
series [^, |62|. We shall start with the so-called regular case, i.e., when 
the series on the r.h.s. of (0) has nonzero radius of convergence. Then, by 
virtue of Theorem]^, our method will sum J2'^=o^^ to 1/(1 — z) in the whole 
complex plane except the ray [l,cx3) (=MLS'[1/(1 — z)]) in contrast to the 
Borel method (see Fig. ^. 

In the next we shall consider moment constant summability method with 
the momemt sequence {/i(^)}5^o, 

jj,{n) := / exp(-expt) t" dt, (14) 

Jo 

i.e., with the moments in the Stieltjes form. Firstly we shall prove that 
the moment constant method solves the above difficulties (A) and (B) of 
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the Borel method, and, in the next, we shall give a whole family of such 
methods. Indeed, one of our main results is that if f{z) is the principal 
branch of an analytic function regular at the origin, where f{z) = J^o-n^^, 
then the integral I{zo), 



POO J^ 

I{zo) = exp(-expt) ^a„ 



>"*'"*. 



n=o f^{n) 
converges if and only if Zo G MLS{f). If Zo G MLS{f) then 

f{Zo) = I{Zo). 

Thus, the method is analytic and regular. The convergence is absolute and 
uniform in any bounded subset of MLS{f) with nonzero distance from the 
boundary of MLS{f), and one can also differentiate inside the sign of inte- 
gration. As for the property (B) the method is analogous to the Lindelof one 
p. In contrast to the Borel method the method does not only see beyond sin- 
gularities and may be used for localization of critical points of some physical 
models. We shall also show that the issues of the Borel summability method 
consist in that it is based on the asymptotic properties of the Mittag-Leffler 
function approaching zero in the complement of some sector with nonzero 
opening angle only [^ |106|| , while the new method (p!4D is based on an entire 



function which approaches zero in all radial directions except for one. In an 
application to the Rayleigh-Schrodinger perturbation theory, theory of linear 
operators in Banach and Hilbert spaces is used [|109|| . 
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3 New results 

3.1 Regular case 

From some general considerations we expect that if a moment constant sum- 
mability method exists solving problems (A) and (B), then the moment 
sequence has to grow like (Inn)" (see Remark 5). However, the problem of 
finding a nondecreasing function x(^)) 

fi{n) = / rdxit), 



for a given moment sequence {/i(^)}5^o (the Stieltjes problem [0, |119|| ) is 



very complicated. Throughout the paper we shall only consider absolute 
continuous functions x(^)- Instead of a moment sequence, we shall pre- 
fer to choose a suitable weight function. Such a weight function could be 
exp(— expt). The reasons are that the moment constant summability me- 
thod with the moment sequence {^(n)}^^^ where /i(ra) are given by (|I^) 
plainly belongs to the moment constant summability methods discussed in 
[|6^ , and as such it is regular. The next reason is yielded by the following 
Lemma. 



Lemma 1 .■ The function 

/■oo 

/i(s) := / exp(— expt)t''(it (15) 

Jo 

1) has a meromorphic extension onto the whole complex plane with simple 
poles on the negative real axis; 

2) its asymptotic behaviour for s tending to infinity, \ args |< vr, is governed 
by the saddle point only and 

/i(s) ~ [2iTw{s)w'{s)]^''^ exp[— expw{s) + s\nw{s)]. (16) 

The saddle point w{s) is given in the complex plane by the equation 

{expt)t = s. (17) 

3) The function fi{s) exhibits no zero for Re s > —1. 
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Before the proof we give some comments on the solutions of (|T^) for s com- 
plex. In this case we shall always assume that the s-complex plane is cut 
along either the real positive either the real negative axis. Correspondingly 
a sheet in the w-complex plane is chosen which is bounded either by the 
real positive axis from below and by a curve + arctan(0/w) = 27r when 
Rew > 0, (f) + arctan(0/w) = vr when Rew < 0, from above, (p being Imw, 
< Imw < 2tt. ; either by the curves (f) + arctan((^/w) = tt when Rew > 0, 

< 0, and + arctan(0/u7) = when Rew < 0, < 0, — vr < Imw < vr. 

Proof : 1) The integral on the r.h.s. of (^) converges absolutely in the 
whole complex halfplane Res > — 1 and defines there an analytic function. 
To perform an analytic continuation we make use of the standard Cauchy 
procedure. Being an entire function, exp(— expt) possesses a power series 
expansion, 

oo 

exp(-expt) = XI cX, 

n=o 

with infinite radius of convergence. Then splitting the integration in (|T5|) 
into two intervals, one obtains 

oo p -i-n+s+l poo 

M^) = E , "° ^n + / exp(-expt)f rft, (18) 

rrl in + s + 1) Jt„ 



n=o 



where the r.h.s. of (jT^) is defined for all s except for s = —n — 1, thereby 
providing an analytic continuation of the integral on the r.h.s. of (p!5D onto 
the whole complex plane after one sets to = 1. 

2) For s tending to infinity and | args |< 7r/2 the asymptotic behaviour 
of fi{s) can be calculated directly using the saddle point technique on the 
integral on the r.h.s. of (|15|). The contribution Vg of the saddle point (|1^ is 
calculated in the standard manner [|104| , |105|| . One rewrites integrand of (p!5|) 



as an exponential function exp/i(t), where 

/i(t) := -e* + slnt, (19) 

and looks for critical points of h{t). This means that one solves the equation 

h'{t) = -e* + s/t = 0. (20) 
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Any solution of ( pOD is called the saddle point and its contribution to (|T5|) is 
given by the formula 



where tg denotes the value of a given saddle point. Under some additional 
conditions one can prove that for s —>■ oo the value of integral (|T3|) is com- 



pletely determined by the contributions of saddle points ||104| |105|| . In our 



case one finds that ( poD has an unique solution in the cut complex plane, 
which we shall denote by w{s)). Since w{s) expw{s) = s one finds that 

^'^') = Ms)f ] ^^iV ^i^Wi^) = -7777^ ■ (22) 

e^'^^>[w{s) + 1) h"[w[s)) 

Thus the contribution of the saddle point w{s) is given by 
Vs = {2nw{s)w'{s)y^^ exp{-e'"^''> + w{s)lns) 



and the asymptotic behaviour of /i(s) is completely governed by it [|104| , |105 
so that 

/i(s) ~ i27iw{s)w\s)Y/^[exp{-e'^^'^ + w{s)\ns)]. 

For s ^ oo and | args |> 7r/2 the situation is quite different, because the 
contribution V^ of the end point of the integral on the r.h.s. of ([18|) ||104|, 105 



Ve = -tlexp{-expto) /s (s -^ oo), 

is far from being negligible and in fact is greater than the contribution Vg 
of the saddle point which approaches zero when s tends to infinity, and 
7r/2 <| args |< vr. However, we can prove that the contribution Vg is nothing 
but the sum on the r.h.s. of (|T^) in the large s limit with opposite sign, so 
that these contributions cancel each other and only the contribution of the 
saddle point Vs survives. In fact, 

oo j.n+s+1 oo j.n 

y Cn— = iC/s) y Cn /° ,, • (23) 

to n + s + 1 ^" ^ 't'o^ + in + l)/s ^ ' 

The sum on the r.h.s. of (|23D converges uniformly in s, <| args |< vr — e, 

oo j.n oo 
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where 6 is some constant whish does not depend on s, so we have 

oo oo 

]hn J2 cJll[^ + (n + l)/s] = ^ c^ = exp(- exp t^), 

n=o n=o 

and 

oo 

E C„tr^+V(« + 5 + 1) ~ -K (5 ^ oo). 

n=o 

To prove 3) suppose that for some So, 

/•oo 

Ai(so) = / exp [—exp (exp t) + (sq + l)t](it = . 

J oo 

However, this imphes that for any 6, 

/•oo 

/ exp{— exp[exp(t + 6)] + (sq + 1)0 dt = , 

J OD 

what means on the other hand that 

/•oo 

/ exp [—exp (exp t) + {so + 2)t][— exp(expt)] dt = 0. 

J oo 

Repeating this procedure m-times one gets that 

/•oo 

I{m) = / exp[-expt+mt+(so+m) \nt]A-"'[A{-l+d/dA)]"'^\-A) dt = 

Jo 

where A = expt. Note that the term D{t), 

D{t) := A-'^[A{-1 + d/dA)Y'''\-A), 

is bounded on t G (0, oo) and nonzero for sufficiently large t. On the other 
side, provided m is sufficiently large we may use an asymptotic formula to 
evaluate I{m). Thus, 

I(m) ~ [2iTu'u/{u + 1)]^ exp[— exp-u + mu + (sq + m) \nu]D{u) = 0, 

since u{m) is the solution of the equation 

-e* + m + {so + m)/t = 0, 
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which behaves hke 

ulm) r^ ln[m + [so + m)/ In m] {m—>-oo). 

O 

If one uses the approximate solution to the saddle point equation (|11 

W7(s) = ln(s/lns) +Inlns/lns + 0[(lnlns/lns)^] (s -^ oo), (24) 

then one finds the following approximate asymptotic behaviour of /i(s), 

/i(s) ~ (27rlns/s)^/2[exp-(s/lns)][ln(s/lns)]' (s ^ oo) , (25) 

Some words about a numerical implementation of the method. Note that 
the method can be implemented numerically as well, since the solution to 
( P!7D is very quickly determined by the following recursive formula: 

w{n) = ln{n / ln{n / ln{n / . . .))), (26) 



and the asymptotic formula ([16D very accurately reproduces the value of the 
integral (p!5| ) (in fact, up to ten orders by a very simple algorithm). One also 
finds that and when s tends to infinity and |args| < tt the relation (p^), 
w'{s) ~ 1/s — l/sw{s), still holds. 

Remark 1 : Note that the proof of Lemma can be directly adapted to the 
case of the weight function exp[— eA;(t)], where ek(t) is the fc-fold exponential 
function, 

efc(t) := exp(exp(...(expt)...), 

k being an arbitrary non-negative integer. The saddle point is now given by 
the equation 

ek{t)ek-i{t)...eh = s. 

In the case oi k = one obtains a simpler proof of the asymptotic properties 



of the gamma function than the one usually known ||106|| . For k > 2 one ob- 
tains another moment constant summability method /i^ solving the problem 
(B) (see below). They give nothing new in the regular case, but they play 
quite an important role in the singular case. Considering the integral on the 
r.h.s. of (|T5|) as the Mellin transform, one could adapt Lemma |1| to study 
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the Mellin transform of some class of entire functions. (} 

Now, our main problem is the determination of asymptotic properties of 
the moment function F{t), defined by 

The reason is that for | z |< 1 

/ exp{-expt)F{tz)dt = ^ z" = 1/(1 - z), (28) 

•^° n=o 

due to the regularity of the /i-method [^]. Because of Lemma |1| the function 
F{t) is an entire function. One expects that 

max In | F{t) |~ exp(constexp | t |) (t ^ oo) 

\t\=const 



105| , |107|| . If F(t) grew no faster than exp \ t \^ when t approaches infinity, 



and I argt |> 0, where A is a constant, then (^) should converge absolutely 
for every z G MLS'[1/(1 — z)], the Mittag-Leffler star of the Cauchy kernel 
1/(1 — z) providing in such way an analytic continuation of I]^o ^" from the 
unit disc onto the whole MLS'[1/(1 — z)]. 

Let f{z) be the principal branch of an analytic function regular at origin. 
Then 

f{z) = (l/2vr0 / -/^ du 
Jc U[l — z/u) 

= (l/27r2) / ^-^ du r exp{-expt)F{tz/u) dt (29) 

Jc U Jo 

for every simple contour C such that no singularity of f{z) lies on C or inside 
it. If in addition the contour C is such that {z/u \ u G C} C MLS'[1/(1 — z)], 
then we can invert the order of the integrations in (pQ]), to obtain 






n=o K^) 
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l-oo °° (tz)^ 

/ exp(-expt) Va„^r^(it, (30) 
Jo frl a{n) 



a„ being the Taylor coefficients of f{z) at origin. It is a short exercise to 
show that the set of z for which the /i-sum ( pOD exists is just MLS{f). Hence 
to show that the /z-method provides an analytic continuation of the Taylor 
series of f{z) at origin onto the whole MLS{f) we have to prove that F{t), as 
defined by (|27|) , does not grow faster than exp(| t |^) when t tends to infinity 
and I arg t | > 0. One can also deal with the case where some finite number of 
singularities of f(z), but no branch point, are elements of Co (interior domain 
with respect to C) provided the condition {z/u \ u G C} C MLS'[1/(1 — z)] 
is satisfied. The result is 

i\z) = -^Res[f{u),Zs]/{zs- z)+ exp(-expt)^a„^-— dt, (31) 

where 

an = a„ + Y.Res[fiu),z,]/z:^\ (32) 

s 

and J2'g in ( pT| , |32D runs over all singularities Zg of f{z), which are elements 
of Co (see also the next section). As an exercise one can justify the relations 
(PT| , P^ for the Cauchy kernel 1/(1 — z), where a„ = so that the integral on 
the r.h.s. of (|3TD vanishes. 

Let us turn now to the study of asymptotic properties of the function 
Fit). 



Lemma 2 Let {/^(n)}^^^ be the Stieltjes moment sequence generated by the 
measure exp(— expt) dt (fZ^. Then the function F{t), defined by (j^^, 

oo j.n 



^0 /^(^) ' 

is an entire function with the following asymptotic behaviour at infinity : 
1) For I Imt \< 7r/2 the asymptotic behaviour is determined by the saddle 
point s = (expt)t, of the Euler-Maclaurin integral representation of F(z), 



F{t)= ds + 0{\t\'') (t->oo) 

Ja jJ[S) 

a being some constant, — 1 < cr < 0. Therefore 

F{t) ~ exp [exp t + t + ln(t + 1)] (t ^ oo) . 
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2) For I Imt \> 'n 12 the asymptotic behaviour of F(z) is still determined by 
the saddle point whenever its contribution prevails the contribution of the end 
point of integration. Otherwise, as well as for \Imt\ > it, 

I F(t) |<0(| t I") (t^oo). 



Proof : By Lemma |T] /i(s) exhibits no zero in the right complex halfplane. 
Then by virtue of the Euler-Maclaurin sum formula ||104] , |105|| we arrive at 
the ensuing integral representation of F{t), 



oo gSlnt 



F{t)= —-ds + 0{\tn (t^oo), (33) 

where a is some noninteger from the interval (—1,0). To find an asymptotic 
behaviour of F{t) at infinity we shall use the standard saddle point method as 
in the proof of Lemma 1. In the present case one does not know an explicite 
form of the function h{s) (|19|) which determines the position of saddle points. 



Nevertheless the position of the saddle point of (|33D can be found exactly 
for sufficiently large t. As one expects and as we shall show in a moment, 
function h{s), and thus equation for the saddle point is determined by the 
asymptotic behaviour of /x(s), provided t is sufficiently large. Let as above 
w = w{s) be the solution of (|I7D , i.e., 

w{s)exp[w(s)] = s. (34) 

Due to Lemma |1| the asymptotic behaviour of fi{s) at infinity is 

/i(s) ~ [2-Kw{s)w'{s)]^''^ exp{— exp[w{s)] + s\nw{s)}. (35) 

This in turn determines an asymptotic form of the function h{s), 

h{s) = exp[w{s)] — slnw{s) + slnt, (36) 

and hence a solution of the equation 

h'{s) = e'^w'is) -\nw- sw'{s)/w{s) + Int = 0. (37) 

By using the relation (|1^ the last equation is essentially reduced to 

h'{s) = - In w + Int = 0, (38) 
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and the saddle point is t = w{s). From (|38| ) one finds that h"{s) = —w'{s)/w{s 
at the saddle point. Finally, from (0) one obtaines the contribution Vs of 
the saddle point, 

where 

u{t) ■=expt + t + \n{t + l). (40) 

If I Imt \< 7r/2 then one can justifies that the asymptotic behaviour of F(t) is 
determined by the contribution of the saddle point only. If 7r/2 < | Imt \ < vr 
then the asymptotic behaviour is still determined by the contribution of the 
saddle point provided it prevails the contribution of the end point of the 
integral (0). For example this happens if one stays on special curves which 
approximates y = tt/2 {y = — vr/2) from above (below) if a; ^ cxo (see also 
below section p.3|) . Otherwise, and for |/r7it| > tt as well, the contribution of 
the end point s = a of the integral ( pB]) prevails and |-F(2;)| tends to zero. If 
I Imt I > 71 the contour of the steepest descent does not exist since the saddle 
point does not lie on the first sheet of the multivalued function w{s). Another 
point of view is that we do not have any saddle point of the integral (^) 
since w{s) in (|35D is only defined for | args |< vr for which | Imw{s) |< tt. 
Hence, the relation ( pHf ) cannot be satisfied for such t. In both cases, however, 
one can show that the asymptotic behaviour of F{t) is determined by the 
end point s = cr of the integral (p3D, i.e., F{t) ~ 0{\ t |'^) for t tending to 
infinity and | Imi |> vr [ 104 , 105|| . For | argt |> this can be examplified 



independently as follows. By virtue of the analytic properties of /i(s) and 
its behaviour at infinity ( pSj) one can rotate the contour of integration from 
the real axis to the ray (cr, a + iqoo) , where q = sign (argt). Let argt > 0. 
Then the integral J, 



-ioo ^slnt 



ds, 



fi{s 
can be majorized in the subsequent manner, 

I / |< e"'""!*! / , , , , \ds\< 2const e"^''!*!, (41) 

Ja I fi{S) I 

where s = a + iy, because on the ray (cr, a + zoo) 

I l//^(-s)l ~ exp [const {y/\ny)] {y ^ oo). (42) 
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Thus I decays at worst like | t l'^ when t tends to infinity and the proof is 
finished. O 

Remark 2 : Note that the position of the saddle point is 

w{s) = t 

also in the case of Stieltjes moment sequence {/ifc(n)}^o generated by exp[efc(t)], 
where k is an arbitrary nonnegative integer. Now, w{s) is the solution of (^) 
with expw; replaced by {d/dw)ek{w) (see Remark 1). The contribution of the 
saddle point is up to the Gaussian integral around the saddle point 

Vs = exp[ek{w)] (t — > oo). 

Analogously as in the above case oi k = 1 , Fk(t) , defined by 

oo j.n 

n=omn) 

where A; > 2, is also polynomially bounded when t tends to infinity and 
I argt |> because of 

I l//"fc(s)| --exp{const[y/\ny...\nk-i{y)\nk{y)]} {y -> oo), 

on the ray (cr, a + iqoo), In^d/) being the fc-fold logarithm. It means that 
such methods also solve the problem (B)! The principal difference between 
the Borel {k = 0) and /ifc-methods with k > 1 consists in the fact that 



for the Borel method the integral on the r.h.s. of ( PD may not converge if 
I argt |> 0. 

The main result of this section is the following Theorem. 



Theorem 2 .■ Let f{z) be the principal branch of an analytic function regular 
at origin, 

oo 

f{z) = Y: anz\ 
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Then the integral I{zo), 

I{zo)= exp(-expt) ^a„ ° dt, (43) 

Jo n=o l^W 

converges if and only if Zo G MLS{f), where fi{n) are defined by (\m). If 
Zo e MLS{f) then 

f{zo) = Hz,). 

The convergence is absolute and uniform in any bounded subset of MLS (f) 
with nonzero distance from the boundary of MLS (f) and we can differentiate 
inside the sign of integration, 

POO °° z"'~^t'^ 

f\z) = / exp(-expt) y^'^^'ri — r~T' ^^■ 



Provided f{z) has no branch point on or inside C, then the relations (|?7|,[^ 
hold. 

Proof : Note that the integrand of (^) is an entire function, so that it is 
well defined. The result follows immediately from Lemma ^ and the relations 



(m®. 



Remarks : MLS{f) is invariant with respect to differentiation, i.e., MLS[f{z)] 
= MLS[f'{z)]. If f{z) is as in Theorem 0, then f'{z) also satisfies the con- 
ditions of this theorem, and we have another representation of f'{z), 



ra-l 



zty 

' dt. 



f>( \ ( +\^ ^^ 
f [z)= / exp(- exp t) 2^ nan—, 

Jo „^i ^^{n - 1 

i.e., the moments ^{n) can be in some sense shifted. <) 



Let f{z) be as in Theorem ^ and regular on the real positive axis, for 
instance. Then Theorem ^ enables us to give an analytic continuation of 
the Taylor series of f{z) at origin on the whole axis. In fact the following 
Corollary holds. 
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Corollary 1 ; Let 

1) f\z) he the principal branch of an analytic function regular at origin; 

2) f{z) he regular on the real positive axis. 
Then 



f{x) = / exp(-expt)^, 

Jo „_Q 



yXt)"" 



/i(n) 

for X G [0, oo). The integral converges ahsolutely and one can differentiate 
inside the sign of integration. 
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3.2 Some new results on the Borel summability me- 
thod 

In this section we wish to prove an analogue of Theorem ^ for the Borel 



method. Some results are exposed in [^ , but a detailed study of the region 
of convergence of the Borel integral is missing. The same role as the function 
F(t) has in the /i-method discussed above plays the generalized Mittag-Leffler 
function in the Borel method. As we have mentioned above the Borel method 
is one member of the moment constant summability methods. Its weight 
function is {l/a)t^^^"^~^ exp{—t^^°'), and the moments of the Borel method 
are in general 



(1/a) / t^'^/")- 



^e"* f' 



dt = T{an + /?), 



where F is the usual gamma function, and a, (3 are positive constants, < 
a < 2 (note that the Borel method with a = (3 = 1 can be viewed as a 
special case Ho of the /i^-methods discussed in Remark 1). The generalized 
Mittag-Leffler function is the function 



EaA^) = E 



T{an + (3) 



[|106|| . We wish to show that the problems (A) and (B) of the Borel me- 
thod are results of the fact that Ea^z) is unbounded in some sector-like 
domain with nonzero opening angle. Indeed, its asymptotic properties when 
z approaches infinity are as follows |6^, |106|| : 
a) 



EaM ~ E -^'Vri/^ - an) if I aig{-z) |< (1 - a/2)7r, 

n=l 



b) 



E^ 



A^) 



;i/«)E4-''e 



if 



axgz \< 7ra/2, 



where tm = ^i/"e^'^*'"/" , and the sum runs over all m such that —na/2 < 
axgz + 2TTm < TTa/2. Because of these asymptotic properties the integral 

POO 

/ t^-'e-'E^A^ndt, (44) 

Jo 
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which is for | z |< 1 equal to 1/(1 — z), does not converge in the whole 
ML5'[1/(1 — z)] but only in some domain Ba[l/{1 — z)], 

i3„[l/(l - z)] := {z I i?e(2)^/" = r^/"[cos{e/a)] < 1}. (45) 

Note that ;Bq,[1/(1 — z)] does not depend on /5. Due to the relations (p^,^) 
such a representation of the Cauchy kernel determines the representation of 
an arbitrary function f{z) regular at origin. Let C{zo) be a contour given by 
the relation 

r = ro{cos[{e - 0o)/a]} where \ 6 - Oo \< na/2. (46) 

Let us draw the contour C{zs) for each singularity Zg of f{z) and let us dis- 
card from the complex plane the domain Co{zs) closed up by the contour, 
for which the sign = in (^6|) is replaced by >. Then the notation Ba{f) is 
adopted for the simply connected region of the complex plane containing the 
origin. By the construction Ba{f) is a starlike region containing, in any case, 
the disc of convergence of the Taylor series of f{z) around the origin. 



Theorem 3 .■ Let f{z) be the principal branch of an analytic function, re- 
gular at origin, 

oo 

f{z) = J2 anZ""- 
n=o 

Then the integral I{zo), 

Jo ^ T[(3 + an) 

converges if and only if Zo G Ba{f)- If Zo G Ba{f) then 

f{Zo) = I{Zo). 



The convergence is absolute and uniform on any bounded subset of Ba{f) 
with nonzero distance from the boundary of Ba{f)- For the first derivative 
f'{z) of f{z) the following representation 



oo n—lj-an 



f'{z)= t^-'e-'Y.nan "' , dt (48) 

Jo frl Lip + an) 



r(/? + an) 
holds. 
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Proof : Proof of Theorem is nothing but a shght modification of the proof 
of Theorem ^ to the case of another Stieltjes moment sequence. The double 
integral on the r.h.s. of ( PD| ) converges if and only if for a given z^ a contour 
C exists such that {zo/u \ u G C} C Ba[\-/{)- — z)]. Interchanching the order 
of integrations in (pOf ) one obtains the relation (|47|) if and only if: 

a) the origin and Zo are elements of Co ; 

b) no singularity Zg of f{z) is an element of Co U C; 

c) Vm G C, where u = pexp{i(j)) and such that | 9o — 4> l< 7ra;/2, the 
inequahty 

Re {zojufl'^ = {ro/pY^" cos[0o - 0)/a] <l-6 (49) 

holds, where 5 is a positive constant. 

The first two conditions are obvious; the last one follows from the fact 
that the integral (H^f) converges if and only if 



Re z^l'^ <\-b. 

One can easily check that a set of all Zo such that a contour C with the above 
properties exists is just the region Ba{f)- In fact, for all z G Baif) such a 
contour exists, li z ^ Ba{f), then inevitably a singularity Zg of f{z) exists 
such that 

r>rs{cos[{9-es)/a]}-'^ when \9-9s\<na/2. 

This, however, contradicts to the properties b) and c) under which 

r, >r{cos[(^-e,)/a]}°>p if \ 6 - 6^ \< 7ia/2. 



The relation ( ^5[ ) then follows immediately from the absolute convergence 
and the fact that under the conditions of the Theorem the integrand of (^) 
is an entire function. O 

Remark 4 : Like MLS{f) the domain Ba{f) is also invariant under differen- 
tiation, Bo,{f') = Ba{f), so that inserting the first derivative f'{z) in place of 
f{z) in Theorem ^ one obtains another integral representation of f'{z), like 
in the previous section. 



oo 



oo / ^j.a\n—l 



lo ^^ T[f5 + a[n-l)] 
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Remark 5 : Note that whenever a < a', then Ba{f) C Ba'{f)- For a ^ 
the domain Ba{f) approaches MLS{f). However, the hmit 

\imT((3 + an) = T((3) 

is trivial (does not depend on n). This fact provides one of arguments for the 
moment sequence of the moment summabihty methods summing the Taylor 
series of f{z) in the whole MLS{f) to grow like (Inn)", because 

/i(n) ~ (Inn)" = o[r(/3 + ara)] (n ^ oo), 

for all positive a and (3. Another way is to consider the weight function 
for the Borel moments which contains a factor exp(— t^'") and its behaviour 
when a — > 0+. <C> 



Theorem |^ shows that in the regular case the Borel summation is not only 
the Laplace transform (cf. ||108|| ). The integral (|47|) in the Laplace form, 

oo j.Qn 

dt. 



I\zo) = (1/zo)'/" / t^-' exp(-t/^y") Y. an 

Jo „ — „ 



T{(3 + an) 



has a different region of convergence but one is an analytic continuation of 
the other. In fact, adopt the notation L{z) {L°{z)) for the simple contour 
(the interior domain with respect to it), where L{z) is parametrized by u : = 
pexp(i0) in the following way (see Fig. ^: 

a) If I ^ — |<| — 0o |< 7ra/2, where z = rexp{i6), then 

p = r{cos[(^-0)/a]}, (50) 

±0o being determined by the equation 

p = r [(cos(^ - (j))/a]'^ = e, 

where e is strictly less than radius of convergence of the Taylor series at the 
origin of the function under consideration; 

b) If I ^ — I > I — 0o I , ±0o being the same as above, then 

p = e. (51) 

The following theorem holds. 
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Theorem 4 ; Let f{z) he the principal branch of an analytic function regular 
at origin. Denote by Da{f) the domain 

DM) ■■= U^eBM)L"{z), (52) 

where L°{z) is defined by (^J5|j. Thus for any z G Da{f), z = rexp{i9), 
there exists z* = r* exp[i6*) G Ba{f) such that z e L°{z*). Let Zq be 

Zo := z exY){—i9*) . 

Then the integral 

l\z):={l/zoY''' t'^-'exp{-t/zl/'')an{e'^*t'')''/T{(3 + an) dt (53) 

Jo 

converges (outside the Borel polygon!) and equals to f{z) . 

Proof : Suppose the corresponding z* G B^if) has been found. Then the 
integral 

POO °o (y*f<^Y 

/ t^'^e-' E a„- ^' > 

Jo „^„ J 



dt 



(1/r)^/" /'°°t^-^exp(-t/r^/°)X]«n A^ ' ^ . dt, 

Jo „_„ 



^o r(/5 + an) 

converges. After the real substitution t — > tr^''" one gets the integral in the 
Laplace form, 

^o^"r(/5 + an) 
which converges Vw G Co (the complex plane of r) such that 

or, what is the same, 

w < r[cos(6'/a)]". 

But the complex plane Co is nothing but the rotated C, i.e., Co = exp{—i6*)C . 
O 



Let us give an example to illustrate the Theorem. Consider the function 
f{z) = 1/(1 — z) and the standard Borel method {a = P = 1). Now, Bi{f) 
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Figure 5: Analytic continuation beyond Re z < 1 . 

is the complex halfplane Rez < 1, and Di{f) = MLS{f). Let us calculate, 
e.g., /(2 + 3i). As the corresponding z* we choose z* = 8i. Hence 6* = it/2 
so we have Zo = 3 — 2i (see Fig. ^. According to ( p3D 



POO j^ 

/(2 + 3t) = (3 - 2ty' / exp[-t/(3 - 2^] ^ 



/•CXD r /I 

,_1 / I . / 1 



{3-2i)-^ exp -t :-i)}dt 



3-2i J } l + 3i 



Theorem ^ shows the principal obstruction for the Borel method to be 
convergent in the whole MLS of function under consideration. It is due to 
the fact that whenever the Borel sum exists for some Zo, then it also exists 
for z e L°{zo) and defines an analytic function there. 

At the end of this section we wish to show that the Borel method can 
represent a function regular on the real positive axis only under assumption 
of its regularity in a larger sector-like domain. 
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Corollary 2 ; Let f{z) be regular in a sector-like domain S^ , S^ := {z 
r > and | ^ |< 7}. Then the integral 



Jo „_„ 



n=o r(/5 + an) 



dt 



converges for z G (0, cxd) if and only if a < 27/7r. If a < 2'y/n, then the 
integral converges on each sector Sn, Sr^ := {z \ r > and \ 9 \< rj}, where 
rj < TTa/2. 

Proof : Only under the condition a < 2'j/a the real positive axis belongs to 

Baif). <> 

Note that all results of the section can be generalized to the case where 
some singularities of f{z), but no branch point, lie in Co- Note also that for 
the class of functions satisfying the hypotheses of the Nevanlinna theorem 
the Borel transform is analytic in a strip including the real positive axis. So, 
applying Corollary |I] one directly finds an expression for the Borel transform 
on the whole axis (in terms of a^*). 

As for the generalization of the strong asymptotic condition one can see 
that by replacing the factor A^! in the above bound of Rn{z) by /ifc(-A) one 
obtains uniqueness theorems for horn-shaped regions, k characterizing the 
sharpness of a horn. The proof is simply performed by combining conformal 
mappings with the Phragmen-Lindelof theorems [|TT| . 
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3.3 Singular case 

As it has been shown above the moment constant summabihty method with 
the moment sequence {/i(^)}^, 

POO 

li{n) := I exp{-exp t)rdt. (54) 

Jo 

solves the shortcoming (B) of the Borel method. In this subsection we shall 
show that this method can be used for the horn-shaped singularity as well. 
We shall only deal with the horn Hji (see Fig. |]), defined as 

Hr := {z I Re uj{l/z) > uj{l/R)}, (55) 

where 

u;{z):=lnF{z), F(^) ;= ^ (56) 

n=o f'K^J 

i.e., roughly speaking, with the region of the asymptotic type (1, 1) [^ ^, 
6|. 



This can be confirmed as follows. By Lemma ^ we know that Reu{z) 
tends to — oo for z ^ oo and Imz > n. On the other hand, within the 
strip Imz < n, uj{z) behaves according to (0) as uj{z) = e^ + z + \ia{z + 
1) provided Reuj{z) does not decrease too fast. Thus within the strip the 
equation Reu{z) = const which defines the separatrix of the asymptotic 
behaviour of F{z) is nothing but 

Reuj{z) = e^ cosy + x + \n\x + iy + 1\ = c, (57) 

c being some real constant. The Eq. (|57D shows that dH^ is symmetric 



under y -^ ±y. Hence we can restrict our considerations to the right upper 
quadrant. In order that Reuj^z) be a constant the exponential term in (0) 
forces y to move very rapidly in the region 7r/2 < \y\ < tt and then to ap- 
proach very close \yo\ = vr/2 from above if a; tends to infinity. One anticipates 
that the boundary dHu of Hr encloses in some sense the strip Imz < 7i/2. 
Indeed, asymptotically, for a; ^ cxd, the boundary of dll]^^ is approximated 
from below and above by the curves 



y±[x) 



- + e "^ [x + \n 



x + l + i(^ + e~''{x±2 ln(x + i; 

(58) 
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Figure 6: The horn Hr (a) and its image Hp} (h) under the mapping z — *■ 

l/z. 

with X G [xo, oo], Xo > 0, since 

Reuj{z)\y=y^(^x) = e^cos|/±(x) + x + In |x + 1 + iy±\ 

, 7rln(x + 1) _„ , , ,^^, 

~ c + a± In — ^-p— ^e ^ (x ^ cx)) (59) 

where a+ = — 1 and a_ = 3. The above considertions confirm our expec- 
tations that the boundary of dHu approximates the origin with zero slope, 
since according to (^,^) , 

y ~ ±2^^ (^ -" ^)' (60) 

where the bared parameters correspond to the mapping z ^ z = l/z. 
From (|57|) it follows that 
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diverges for y — >• as expected. The point Xo at which y{xo) = satisfies 
approximately the equation 

(xo + l)e^° = e" (62) 

and the point Xm at which y{x) takes on its maximum is given by the equa- 
tion, 

[(x^ + l)2 + y^]VV™=e'^+\ (63) 



Both (|62D and (|63D are vahd asymptotically for c sufficiently large. Thus, at 



this limit, 



^o + 1 gx„-a.„ _ gC-c-1 _ g-1^ (g4) 



i.e., as Xo ^ oo(<;=^ c ^ oo) then Xm ^ Xo + I. Therefore, by the triangle 
inequahty. 

If ... dz/z\ < 271 max(. . .) + / I . . . \dr/r, (65) 

JoHr Jo 

where . . . stands for any continuous function bounded on the strip \Imy\ < vr 
and X G [xo, C)o), Xo being a positive real number. 

Transition to more general horns of the type [1,1]) is simply accomplished 
by mapping z — >■ z/rj. The theorems could be also modified for a region of 
the asymptotic type {k,r)) with k > 1. The situation here is, however, more 
involved. The main result of this subsection is the Nevanlinna-like theorem 
for the horn Hji as follows. 



Theorem 5 .■ Let f{z) be analytic in the horn-shaped region H^ := {Re uj{l/ z) > 
ci;(l/i?)}, continuous up to the boundary, and satisfy there the estimates 

N-l 

f{z) = J2 dkz'' + Rn{z) 

k=o 

with 

\ Rn{z) \< A^i{N) \ z \^ (66) 

uniformly in N and z G Hr. 
Then 

oo j.n 
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converges for | t |< 1, and has an analytic continuation to the striplike region 
Si = {t \ dist{t,R+) < 1}, satisfying the bound 

I M{t) \< is:exp[exp(| t \ / R)] (68) 

uniformly in every S^ with k > 1. The analytic continuation of M{t) for 
t G (1, cxd) is given as follows, 

M{t) = ^i F{t/z)f{z)dz/z . (69) 

Furthermore f can he represented by the absolutely convergent integral 

f{x) = / exp(-expt)M(tx) dt , (70) 

Jo 

for any x G (0, R). 



The proof of the theorem is rather comphcated. Therefore, it is divided 
into several lemmas. However, an importance of F{t) is already seen from 
the subsequent argument. If a function f{z) is analytic in the horn Hpt and 
continuous in Hr, then for any x G Hr fl R, 

fix) = -—f f{z)/{z-x)dz 

= TT- f fiz)dz/z [ exp{-expt)F{tx/z)dt (71) 

2m JdHji Jo 

Note that such a representation of f{z) is impossible by the Borel method. 
Note also that unlike the disc Cr in the Nevanlinna theorem the horn H^ 
is not a star-like region anymore. This is in general the main difference 
between the regions of the asymptotic type (0,77) and {k,ri) with k > 1. 
This difference means that with moment constant summability methods one 
cannot recover f{z) from its asymptotic series in the whole horn Hr but only 
for z G Hr n R. Physically this is not, however, a problem, since we are 
expanding in real parameters (couplings). 

Henceforth we shall follow Sokal's strategy of the proof of the Nevanlinna 
theorem |^S[. Lemma ^ provides us with an integral representation of the 
monomials t"//i(n) for any positive n and t > 1, 

1 r t" 

i Fitlz)z''dzlz = — -, (72) 

2711 JdHn fi{n) ' ^ ^ 
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where the integral is taken counterclockwise along the boundary of Hr. The 
formula will be used to express M{t) (an analogue of B{t)) in terms of f{z). 
To find the minimal domain of analyticity of M{t) we shall need a bound 
on F^'^\zo) for n tending to infinity (Lemma |^). Finally, after the Lemma 
^ we shall give all the lemmas together and complete the proof of Theorem ^. 



Lemma 3 ; For any z^ such that Re Zo > and n > 0, 

|i^^"H^o)| < const {n\/i2{n)) 

"^P I^^TTT f 1 + -^) - ^(^) + ^o - 1 + 0{w{n)w'{n))] , (73) 
[w[n) + 1 V w{n) J J 

where Xo = Re Zo- 
Proof : Firstly, 

F^''\zo)/n\ = {l/27ci) j> F{zo + z)/z'^+^dz 

= (l/27ri) j> e:!i.Y>[uo{z + Zo) — nlnz]dz, (74) 

where C is a simple contour enclosing the origin. Due to the properties of 
u{z) (see ( |iO| ) and Lemma |^) if \ImZo\ > tc and if one takes C to be a 
contour of radius 1 then one can show that \F^'^\zo)\ — ^ 0. Hence without 
any restriction we can confine ourselves to the region where \ImZo\ < vr. In 
this region we shall evaluate the integral on the r.h.s of (^) by the saddle 
point technique. The saddle point of the integral is a critical point of the 
function 

h{z) := 6"^°+"^ + Zo + z + \n{zo + z + l)-n\nz, (75) 

i.e., a solution to the equation 

h'{z) = e"°+^ + 1 + l/{zo + z + l)-n/z = 0. (76) 

By direct comparison of (|7^) with the defining equation for w{s) (^7^, w{s)e'^^'^^ 
= s, one anticipates that for n sufficiently large the solution v{n) of ( [76|) will 
be "very close" to w{s). Therefore we shall look for the solution v{n) of (ff^) 
in the form 

v{n) = w{n) — Zo — 5{n), (77) 
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where 6{n) is an unknown function to be determined. Such parametrization 
will be shown to be justified as at the end of our calculation we shall obtain an 
asymptotic expansion of S{n) for n -^ oo, according to which the dominant 
term —Zo/w{n) ~ — Zo/lnn -^ 0. To prove this we shall slightly rewrite 
Eq-dH), 

v{n) ( e^°+-W + 1 + — — i ) = n, (78) 

V v{n) + Zo + lJ ^ ' 

and take logarithms of both its sides. Upon the substitution v{n) = w{n) — 
Zq — 5 and some manipulations (after one has expanded small terms in loga- 
rithms) one gets, 

Inw V w — 5 ^ 1 — ~ In n. (79) 

WW e^ e^{w — Zo — 0) 

By using the defining relation (p!7| ) for w{s), \nw{n) + w{n) = n, and sub- 
sequent multiplication of both sides of the last Eq. by we~^ /{w + 1) one 
obtaines, 



w 



S + — ^ e--^ ~ 1 + - + ^^^ ^^^ ■ (80) 



w + lj \ w w^ Jw + 1 

By using ( p2| ) one finds that the last fraction on the r.h.s. of (pO]) is nothing 
but w{n)w'{n). To get 5{n) in an explicite form one again takes logarithm 
of both sides of (IST 



-5 + \n(5+ -^ ] ~ ln(^t^') + "^ ^y^ ^ ■ (81) 



Zo \ , , ,. W + Zo + 5 

—— ~ \n{ww + 

W + IJ w^ 

Dominant terms for n ^ oo in (^) are logarithms. For a moment we shall 
parametrize 6 + Zo/{w + 1) = ww' + e{n)ww', where e(n) is assumed to be 
a small number for sufficiently large n. Under this assumption one gets by 
expanding small terms in logarithm on the l.h.s. of (|8l|) that 



w[n) w[n) + 1 

Thus, the assumption that e(n) is small for n ^ oo is justified and one arrives 
at the following expression for 6{n), 



>\2] 



S{n) = \ + w{n)w\n) ( 1 + -^ - /." J + 0[{ww 

[n -^ oo). (83) 
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This is one of the most important relations to prove Theorem |^. The asym- 
ptotic expansion (^) justifies the statement made at the beginning of our 
calculations that S{n) = (9(1/ Inn) -^ provided n tends to infinity. 

Finally, by using (|83|), one arrives at the following expression for the 
solution v{n) of (^) in the leading order in n, 



v{n) = w{n) — Zo + 



Zo 



w{n) + 1 



+ w{n)w' {n) + 0[w' {n)] 



{n -^ cx)). (84) 



Note that z^ + v{n) approaches the real positive axis when n -^ oo , and 
a part of the contour of the steepest descent nearby the saddle point is 
approximately a segment of the circle centered at z^. After the saddle point 
evaluation of the integral one finds that 

F^'^\zo)ln\ ~ [i;(ra)f'(n)/27r]^/^exp[u;(zo + t;(n)) -nlnw(ra)]. 

Now we can calculate the behaviour of n(ln w{n) — In v{ji)) and ^.^o^^^"^) — 
^w{n) fQj, ^ _^ J5Q_ Both expressions are needed to establish the bound (^) 
o^F^'^\zo). 



w 



nC[nw(n) —\nv(n)) = nln— = nln 

V 



fo+v{n) ^ I ^ l/v{n) 

pw{n) 



nln e-^ + e""' 



w — Zo — 6 



n[-(5(n) + C(e-'"^"0] =n 



~ n 



Zo 



{e-' - 1) + e-'" + 



win) + 1 



w — Zo — 5{n) 
w{n)w'{n) — 0{w'{n)) 



(85) 



^Zo+v(n) _ ^w{n) 



e«'(n)(e-'5 



n 



l) = e 



w(n) 



win) + 1 



w{n)w'{n) — 0{w') 



Win) 



w{n) + 1 

Thus, when we use that e""*^"^ = n/w{n) 

F^^^\zo)li{n) 



— w{n)w'{n) — 0{w'{n)) 



ivv'ww'Y'^ 



n! 
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exp 



ff,v+zo _^ ^ _^ 2:0 + \n{v + Zo + I) - n\nv - e"" + n\nw\ 



exp < , ° _, 1 H -— - w{n) - nw{n)w'{n) 1 H -— + -Zo + C(w(n)w'(n)) 

[w(nj + l y w[n) J \ w[n) J 

{n^oo). (87) 

The factor e"""*^""^ in (^) is produced by (ff'ww')^/^ as one can check by 
using (H) and (H). O 



Lemma 4 ; Lei x & R and z & C he independent variables and let Jx{z) be 
the integral 



y^^ exp{-expt)rdt. 

^„ uin] Jo 



Jx{z) := / exp[— exp(t/x)]F(t2;)(it 

Jo 

to Kn) 
Then we have 

(p Jx(l/z)z'^dz/z = —-— exp(-expt/x)t'^dt, 

2-Kl JdHa fl[n) Jo 

where the integral is taken counterclockwise along the boundary dH^ of H^. 

Proof : The essence of the proof is to show that the integration (^) of 
F{tz) does not spoil the asymptotic behaviour of F{z) too much. Note 
that x" /q ^exp(— expt)t"'(it behaves hke exp[— exp(l/x)]/(nx) provided n 
is sufficiently large (the asymptotic is given by the end point of integration 
here). Thus one anticipates that x will only produce an overall factor to 
the asymptotic behaviour of Jx{z) when z tends to infinity. This can be 
confirmed as in Lemma |^. One uses the Euler-Maclaurin sum formula and 
finds that for | Imz \< vr the asymptotic behaviour is governed by a saddle 
point (provided one stays on an appropriate curve for 7r/2 <| Imz \< vr (see 
the proof of Lemma |^)), which is here determined by the equation 

Inw(s) + 1/s = Inz. 
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Thus, the contribution Vs of the saddle point is 

Vs = exp{exp z[exp{~ exp{—z)/z)] + 1} (^ ^^ C)o). 

X 

By analogy with Lemma |^ one proves that for z tending to infinity and 
I Imz \> tt/2 + 6, 6 0, function Jx{z) tends to zero like F{z) does. The 
proof of the statement of the lemma is then trivial. It amounts to using 
the Cauchy integral formula. Note that the integral from Jx{z)/z along a 
segment of dH]^^, which starts at infinity and terminates at some Zo on the 
contour, converges even absolutely. O 



3.3.1 Proof of Theorem 5 

i) Under the hypotheses of Theorem |^ one easily proves that the series (^) 
converges for | t |< 1. Let us consider the integral d{t), 

d{t) ■■=^<f F{t/z)fiz)dz/z, (89) 

where t > L Two remarks are in order. In contrast to the Nevanlinna 
theorem one cannot use the integral on the r.h.s. of ( pQ]) for < t < 1, as it 
is not possible to satisfy both conditions that the contour of integration in 
( |72D be the contour which tends to zero on the boundary of tH^ for some r, 
and at the same time lie in Hr- Whenever t ^ R the integral is identically 
zero (by virtue of Lemma 0). Hence, it cannot yield an analytic continuation 
of d{t) for t ^ R. 

From the properties of F{z), it is immediately seen that the integral on 
the r.h.s. of ( pQ]) converges for t > 1 absolutely and uniformly on any closed 
subset of (l,oo). For t = 1 the integral converges by virtue of the Abel- 
Dirichlet lemma. Hence, d{t) is a C°° function on the interval (1, oo) and 
possesses the right derivatives at the point t = 1. 

To prove that the series (^) converges at t = 1, we make use of Lemma 
^ and rewrite d{l) as follows, 

N-l 1 

^(1) = E ^ + ^ / F(l/z)R^{z)dz/z. (90) 

^ /i(A;) 27rz JdHR 
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As the integration on the r.h.s. of ( pO] ) runs along the boundary of Hpi on 
which Re uj{l/z) = const (= u^l/R)) due to (|65|) the integral can be bounded 
from above as follows, 

— / Fn/z)RN(z) dz/z < (v4/7r)e'^(^/^V(^) / exp(iVlnr) \dz\lr 

2ni JdHn JdH+ 

< 2A e"(i/^) /i(A^) exp(-Ar ln(l/i?)) (l + ^) • (91) 

Here dH^ means that we integrate along dH^ in the first quadrant. Op- 
timalization of the bound on the r.h.s. of (|^) then amounts to finding a 
solution R = R{N) of the equation, 

h{R) :=-NR + exp{l/R) + 1 + — - = 0. (92) 

1 + 1/ R 

If one compares ( p^ ) with (^) then one finds that ( p^ ) is nothing but (^) 
with Zo = 0. Therefore the bound on the r.h.s. is optimalized by the choice 
R = l/v{N). After distorting the contour of integration up dHji with R = 
l/v{N) and using (^,|^) we have. 



-^ / Fil/z)RN{z)dz/z 



< 



< 2Aexp{e'' + l + ln(t; + l) - iVlni; - e"' + A^lnw - U7/2} 

= 2Aexp{-Nw{N)w'{N)-w{N)/2-0{Nw'{N))} 

r^2Aexp{-3w{N)/2) — ^0 {N ^ oo), (93) 

where we have used that Nw'{N) — > 1 if A^ -^ oo. The term —w{N)/2 has its 
origin in the Gaussian prefector (27rww'Y^'^ of the asymptotic of /i(A^). Thus, 
M(l) = d{l). The same is also true for derivatives, i.e., M('')(1_) = c/(")(l+) 
for any n & N. Indeed, if f{z) satisfies SAC then the same SAC will also 
satisfy its derivatives f^"''{z) (may be in a horn Hr with different R). So 

dWfU) = ^ c,//i(fc) + —i Fil/z)R^;;\z)dz/z, 
, „ Zm JdHn 



k=o 



where Ck = {k + n){k + n — 1) ...{k + l)ak+n, and the integral can be estimated 
in the same manner as above. 
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To determine the ininiinal region of analyticity we express d^^^ (t) as fol- 
lows: 

rf(")(t) = n\a„/fx(n) + — { F^''\t/ z)Rn+i{z)z-''dzl z. 

27r« JdHn 

Now, using Lemma |], one finds that 

I rf(")(t) |<constF(t/i?) n\ 

e^v[{t/R){n/\nn + n/\n^n) + 2\Yiw{n) + 0{l/w{n))\. (94) 

One may justify that 

OO CX) 

Y: d^^\u){t - ti)7n! = Y. d^''\t2)ih - t)7n!, 

n=o n=o 

whenever 1 < ti < t < t2, and | t — tj |< 1, where i = 1,2. 
Thus M(t), 

is shown to be analytic at least in the striplike region Si (see Fig. |^). Bound 
( |68D on M{t) follows immediately from the relation (0). 



ii) To prove the relation (70), let us consider the integral 

1 r°° 

- / exp[-exp(t/x)]M(t)rft, (95) 

X Jo 

for any x G (0, R). The integral, 

/ exp[— exp(t/s)] f f(z)F(t/z) dz/z dt, 

Jl 2'Ki JdHn 

is absolutely convergent. Therefore, the r.h.s. of (|70D (the integral (p5D) can 
be recast into the form, 

/(^) - 77^ / Ullz)f{z) dz/z 

+ - y -P- / exp[-exp(t/x)] r dt. (96) 
X f^o Kn) Jo 
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Figure 7: Minimal region of analyticity of M{t). 

Now we shall show that the last two terms cancel each other. Indeed, using 
Lemma ^, one finds 

— i Ul/z)f{z)dz/z =J2-^ exp[-exp(t/x)] T dt 
2m JdHR ;^ yu(n) Jo 

+7r- / MI/z)Rn{z) dz/z . 

Im JdHn 

To show that the last integral vanishes in the limit iV — >■ cx) one optimalizes 
the bound on the integral. To do this one chooses similarly like in, (|6[]93|), 
the contour of integration to be c}ifi/„(7v)- Then, in virtue of ( ]T7| , |2^J3^ , 

— - f J^{l/z)RN{z)dz/z 

27X1 JdHR 

o 



<— e-3-W/2_^0 (iV^O). 



X 



Corollary 3 .■ Let 
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1) M{t) := I]^o ant"'//i(n) has nonzero radius of convergence and can be 
analytically extended on the real positive axis; 

2) I M{t) \< const exp(exp(t/i?)) on the axis. 

Then 

1 r°^ 
f(x) := - / exp{- exp{t/x))M{t)dt 



X 

converges absolutely for x G (0, R) and 



fix) r^Y.^^nx'' (x^0+). 



Proof : Use the Laplace method to evaluate the integral. <) 
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3.4 Comparison of the new summability method with 
the Borel method 

Using a recent generalization of the Borel summability method [Q and the 



Nevanlinna theorem |jTO[ one can easily find series which violate SAC in a 
horn-shaped region H^ but they are Borel summable (see Example 1 be- 
low). On the other hand one can find series which are formally Borel non- 
summable but satisfy SAC in a suitable kidney-shaped region (see Example 
2). This exemplifies the statement in the Introduction that not the Borel 
(non) summability but the validity (violation) of SAC prevents (indicates) 
the presence of nonperturbative effects which may in some cases destabilize 
a ground state. The first shows a series which violates SAC if summed up 



to a function analytic in the horn Hpi [10|, while it satisfies them if summed 



up to a function regular in a disc tangent to imaginary axis at the origin. 
The main method for constructing them is to use the new AMCSM based on 
the momenta (^4|) as a method for localization of singularities in the com- 
plex plane. Indeed, this method has a nice and rather feasible advantage 
as compared with the Borel one - it does not only see beyond singularities. 
The problem of finding singularities of a function f{z) which is defined by a 
convergent power series is then reduced to the problem of finding the asymp- 
totic behaviour of its moment transform. To find the asymptotic behaviour of 
the moment transform the Euler-Maclaurin formula will be frequently used. 
It reduces the above problem to the problem of finding the asymptotic be- 



haviour of an integral. By our general results |]^, |T0| we know that in general 
the /i-transform may have dangerous asymptotic behaviour (i.e., which may 
prevent summability of its Taylor expansion) at worst on several (possibly 
infinite) radial rays. This should be contrasted with the Borel method - any 
singularity of f{z) results in that the Borel transform of its Taylor series 
grows up very strongly (such that the Borel integral diverges) not only on 
the radial ray which passes through the singularity but also in the whole 
sector bisected by this ray and with the opening angle vr. 



Example 1 Let Si be a power series Y^"^ anQ^ with an = (— l)"/i^(ri), Then 
Si is Borel summable and satisfies SAC in the disc K{0,R) but is not /i- 
summable and violates SAC in the horn H^. 
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Proof: The Borel transform of 5*1 is an entire function and hence it is obvi- 
ously analytic on a strip Sp := {t \ dist(t, i?+) < 1/p} for any p, which is 
the first hypothesis of the Nevanlinna theorem |B^, ^, |TD[. Now we use the 
Euler-Maclaurin sum formula ||105|| , 



B{t) = f e^^'^^^-'^^^fl^rf^ + 0(1 r I) {t -. oo), 

where a is a real number, — 1 < cr < 0. One finds that a saddle point of 
the Euler-Maclaurin integral is at s/w'^{s) = —t. Thus the asymptotic of the 
integral is dominated by a contribution of the end point and we arrive at the 
result 

B{t)^0{\tn (t^oo), 

in any sector with nonzero opening angle | argt |< 7r/2. Therefore the 
hypotheses of the Nevanlinna theorem are satisfied, 5*1 is Borel summable, 
and its Borel sum satisfies SAC in the disc -/^(O, R). 

On the other hand, the generalized yU-transform M{t), 

oo j.n oo 

is undefined since it has zero radius of convergence, i.e., the summability 
method for the horn-shaped region Hpt in this case fails. <) 



Example 2 Let us consider the series S2 = J2'^ a-nO^ with 

a„ = (-lfr(an + l), 

where a = 1 + e, 1>£:>0, and T is the usual gamma function. 

Then S2 cannot be summed up to a function maximal region the analyt- 
icity of which is a disc tangent to the imaginary axis at the origin. Hence it 
is not Borel summable and violates SAC in this region. Nevertheless, it can 
be summed up to a function whose minimal region of analyticity contains a 
kidney shaped region K{X, R), 

K{X, R) = K\0, R) := {z \ z^l^ G ir(0, i?)}, 

having an opening angle = Att with A > a, i^(0, R) being the disc Re{\/ z) > 
1/R, i.e., 5*2 is summable and satisfies SAC in this region. 
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For the proof see |T8|, O . (} 



Eventually, we shall give an example just opposite to Example 1. 

Example 3 Let S^ he a power series with coefficients an = a"'fi{n), where 
a =\ a \ e^^ is a complex number with < 6 < 7i/2. Then the series S3 is 
fi-summable but not Borel summable. 

Proof : The Borel transform of 6*3 is an entire function. The position of a 
saddle point of its Euler-Maclaurin integral representation is given by equa- 
tion 

w = \n{aue'^^+'^^), 

where t = ue^'^, which means that the contribution Vg of the saddle point is 

Vg ~ exp[au exp{i{(f) + 6)){ln{au) — 1 + i{(j) + 6))] {u -^ 00). 

Hence, the Borel integral diverges whenever 6 < 7r/2. On the other hand 
yti-transform M{t) can be calculated exactly, 

°° a 1 

M{t) = Y, ^e 



^oK^) I -at' 



and 5*3 is transparently /i-summable (see |^ |TU 
If 37r/2 > 6 > 7r/2 then one can show that 

Jo n=o ^' "'" n=o f'W 

since the real positive axis belongs to the Borel polygon of M{z). (} 



The above examples demonstrate that there exist series which being Borel 
summable are not /x-summable and vice versa. They also show that /x-method 
[P, |T^, in combination with the Euler-Maclaurin sum formula, is a powerfuU 



tool for looking for singularities of analytic continuation of convergent power 
series in the complex plane, since in contrast to the Borel method /i-method 
does not only see beyond singularities. 
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3.5 Derivation of strong asymptotic conditions for va- 
riety of regions 

To generalize SAC let us start with some definitions. 

Definition 1 Let S be the halfplane Rez > 0, and Cr a disc centered at 
the origin with radius R. Denote S := S\Cr which will be called the base 
region. A domain D of the complex plane C will be called of asymptotic type 
{k,ri), k > 1, if there exist R > efc_i(l) such that D = l/[ri\nk{S)], where 
Ck is the k-fold exponential, eo{x) := x, efc(x) := exp(efc_i(x)), and In^ is 
the k-fold logarithm, lni(2;) := Inz, Inkiz) := \n{lnk-i{z)) . A domain D of 
the complex plane will be called of asymptotic type {0,ri) if D = 1/S^, where 
Sv ■= {z I ^V^ G S}. <> 



To establish SAC for these regions we shall start with the regions of asymp- 
totic type {l,f]). The main result of this subsection is the following theorem. 



Theorem 6 Let f(z) 

1) be regular in a region D of the asymptotic type (1,"^) and continuous 
in its closure; 

2) have in D an asymptotic expansion 

N-l 

f{z) = Y.anz'^ + RNiz). (97) 

n=o 

If R]\f{z) satisfies the bound 

I Rn{z) \< Ap^ I z |^exp[-expw(A^) + A^lnw(A^)] (98) 

with p < rj, uniformly in z & D and N , w{s) being the solution ofw{s) exp[w{s)] 
= s j^, [7^/ , then the asymptotic series ^^ ) determines f{z) uniquely. The 
condition j^OBj ) is strong in the following sense: if it is known that (jP^ is 
only fulfilled with some p > rj, uniformly in z & D and N, then there exists 
a nonzero function with the trivial asymptotic expansion and the asymptotic 
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eries l ^l{ ) does not determine f(z) uniquely. O 



Note that exp{- exp w{N)+N In w{N)} is, up to the factor [2nw{N)w'{N)Y^'^ 
which arises from the Gaussian integration around a saddle point, precisely 
the asymptotic for large A^ of the moment function yu(A^) (see (|1^ above). 
For regions of, roughly speaking, asymptotic type (1, rj) SAC have been estab- 
lished in [^ in connection with a summability method for the horn-shaped 
region Hn. Unfortunately, generalization of this proof to more sharper horns 
seems to be cumbersome. A rather brief derivation of SAC based on a slight 
modification of the Phragmen-Lindelof theorem and which can be straight- 
forward generalized to sharper horn-shaped regions has been given in |]1T| 
and is briefly repeated in Appendix. One has only to replace w{s) by Wk{s), 
Wk{s) now being the solution of 

w{s) exp{w{s) + ... + ek-2[w{s)] + ek-i[w{s)]} = s 

(see 0). One can also consider other types of asymptotic regions which may 
be obtained, e.g., by more involved combinations of scaling z -^ z/rj and 
conformal mappings like z — > 2;^/" of S, as well as some other base regions, 
etc. In general, however, to any given asymptotic region such a moment 
sequence flk{n) will correspond that fulfilling a bound like (^) will ensure 
SAC in the region. We used the above definition of the asymptotic types 
since the maximal region of analyticity of four dimensional renormalizable 
massless field theories has been suggested to be at best just a region of the 
asymptotic type (1, 7), where 7 depends on the first two coefficients /3i, P2 of 
the /3- function and on the definition of the coupling^ ||2^, ^ |9^, ^. This 
means that probably there is no reason to look at the Borel transform of the 
four dimensional renormalizable field theories since the convergence of the 
Borel integral contradicts the horn-shaped region of analyticity [^, |TU|, ^ 



any way that the Borel transform on the real positive axis is defined (e.g., in 



a distributional sense [ p.l(J|| ). In connection with Theorem ^ this means that 
these theories (without UV cutoff) cannot satisfy SAC since their analyticity 
in the complex coupling constant plane is not compatible with the divergence 
of order {n\y of their perturbation theory no matter how small e > is. 



^For 't Hooft's coupling, 7 is just | f3i 
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Now we want to show an important property of the class /C of function 
which obey the hypotheses of Theorem ^ and p. Using asymptotic behaviour 
of Wk{s) [y] and the method of [^ one finds that the SAC we are given 
preserve nonhnear perturbation conditions such as unitarity of the Feynman 
series. Indeed, the following Lemma holds. 



Lemma 5 .■ The class of function K, is closed under product, i.e., if fi{z) 
and f2{z) are two functions from K, then g{z) := fi{z).f2{z) G /C . 

Proof : To prove the lemma it is sufficient to prove that 

N 



lim V fi{N - k)n{k)/n{N) < cx). (99) 

k=o 

From the Lemma |I] one can derive that if N is sufficiently large then for 
q < N/2, 

fi{N - q)/lj{N) ~ exp[-glnln A^ - q/N - {q + q'^)/{N\nN)]. 

Thus, provided one choses a fixed j such that the asymptotic formula to 
estimate /i(s) for s > j can be used, 

N [f+1] 

J2f^{N-k)f^ik)/fxiN)<2 Y: li{N-k)ii{k)/ii{N) 

k=o k=o 

< const ^ exp < — In In A^ 



k>j 



Ink ( Inln/c 

+ A; 1- 



In In A^ V In In A^ , 



The last bound gives transparently convergent series in the limit N -^ oo, 
since its last term behaves like l/[A^exp(A^ln2/ln A^)] . O 

The next two lemmas deal with the important case of perturbation series 
with equal sign and alternating sign of coefficients ||T5|, ^ |2^, ^ |^. To 
formulate the lemmas let us firstly consider a general analyticity region D 
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[n- 


■^ ooj; 


{n- 


-> oo); 


{n- 


-> oo). 



to which corresponds momenta ^k{n). Then one can meet the subsequent 
situations: 

a)\/t > : anf^/flkin) -^ oo 
b)\/t > : a„r//ifc(n) -^ 
c)3t^0 : a„r//ifc(n) ^ ir ^ 



Lemma 6 Let S be a divergent power series with equal sign coefficients. 
Then S violates SAC in any analyticity region D. 

Proof : In case a) the moment constant transform Mk{t), 

oo 

Mfc(t):=^aX//ifc(n), (100) 

n=o 

does not exist. In case b) the moment constant transform is an entire function 
and hence defined for all complex t. It is clear that the maximum modulus 
of Mk{t) for I t |< a; is just Mk{x). From simple relation between the Taylor 
series coefficients and the maximum modulus growth of entire functions based 



on the Cauchy integral formula [p.05|| one can prove that the generalized 



moment constant sum does not exist, since the integral 

/ exp{-ek{t))MkiztP)dt (101) 

Jo 

diverges for all real 2; > 0. The same is also true in case c) since Mk{t) is 
singular on the real positive axis. O 



Lemma 7 Let S be a divergent power series with alternating sign regular 
coefficients a„, i.e., there exist an analytic function a(s) in the complex half- 
plane Res > a such that a„ = (— l)"a(?2). Let a{s) have parametrization 
a{s) = exp(sln6(s)) with sb'{s)/b{s) ~ 0(1) (or o(l)) when s — > 00. Let D 
be the analyticity domain to which the momenta fikin) correspond. Then in 
cases b) and c) S satisfies SAC in D. 

Proof: In case b) the moment constant transform Mk{z) is an entire func- 
tion so the requirement of analyticity of Mk{z) to sum S is satisfied. The 
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Stieltjes moment function fik{s) can be shown to be nonzero in the com- 
plex halfplane Res > a' for a' sufficiently large 0. So, we can again use 
the Euler-Maclaurin sum formula to find asymptotic behaviour of Mk{z) for 
z -^ oo. A contribution of the saddle point is proportional to 

exp{ek[wkis{z))]} ^ exp[efc(-2;)] ~ 0(1) {z -^ oo), 

where s{z) is determined by the saddle point equation Wk{s) exp[— a'(s)/a(s)] = 
—z := — re**^. However, the saddle point does not lie on the principal sheet 
and the contour of integration cannot probably be deformed in such a way 
that the asymptotic behaviour of the Euler-Maclaurin integral be governed 
by it. Instead of that the contour of integration can be deformed in such 
a way that the asymptotic behaviour of the integral can be shown to be 



governed by the end point of integration [|105|] . In any case, however 



Mkiz) ~ 0(1 z n {z ^ oo), (102) 

where a = max{cr, a'}, i.e., the generalized moment constant sum (|101|) ex- 
ists. In case c) Mk{z) is a meromorphic function regular on the real positive 
axis. To see this one may apply //-method f^, |lOl on the series (|10CI|) . The 



/i-transform of Mk{z) is an entire function with alternating sign regular co- 
efficients and this case can be treated as case b). Using the Euler-Maclaurin 
integral representation of the /i-transform Mk{z) one finds that /i-sum con- 
verges Vt > even in a sector and is absolutely bounded therein. Therefore 
/ifc-sum also exists in this sector and satisfies SAC in D. O 
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4 Applications 

4.1 Summability methods and the Rayleigh-Schrodinger 
perturbation theory 

4.1.1 Preliminaries 

As some example of application of the summability methods we shall consider 
the Rayleigh-Schrodinger perturbation theory ||108| , |109|| . We shall restrict 



ourselves to the case of relatively bounded perturbations [|109|| . Let T and A 
be operators with the same domain space X (but not necessarily with the 
same range space) such that D(T) C D(A) and 

Pm|| <a||'u|| + 6||TM||, MeD(T), (103) 

where a and b are nonnegative constants. Then we shall say that A is rela- 
tively bounded with respect to T or simply T — bounded. The smallest lower 
bound bo of all possible constants b in ( |103|) will be called the relative bound 
of A with respect to T or simply the T — bound of A. If b is chosen very 
close to bo, the other constant a will in general have to be chosen very large; 
thus it is in general impossible to set b = bo in ( p.03|) . Obviously a bounded 
operator A is T-bounded for any T with D(T) C D(yl), with T-bound equal 
to zero. 

To consider relatively bounded perturbations is very natural since closed- 
ness, bounded invertibility, selfadjointness, as well as some other properties 
are stabile under relatively bounded perturbation [|109|] . 



Lemma 8 Let T and A be operators from X to Y, and let A be T-bounded 
with T-bound smaller than 1. Then S:=T-hA is closable if and only if T 
is closable; in this case the closures of T and S have the same domain. In 
particular S is closed if and only if T is. 



Proof is rather simple. We shall use the treatment exposed, e.g., in ||109|| . In 
the inequality (|103|) we may assume that 6 < 1. Hence 

- a||M|| + (1 - b)\\Tu\\ < \\Su\\ < a\\u\\ + (1 + &)||Tm||, (104) 

for D (T) . Let us recall that T is closable if and only if u„ G D (T) , m„ —>■ and 
Tun — > V imply v = 0. Applying the second equality of ( |104| ) to u replaced 
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by Un — Um, we see that a T convergent sequence {«„} (that is a convergent 
sequence {«„} for which Tm„ is also convergent) is also S"- convergent. Simi- 
larly from the first inequality an S'-convergent sequence {un} is T-convergent. 
If {un} is S'-convergent to 0, it is T-convergent to so that Tun — ;> if T is 
closable; then it follows from the second inequality of ( |104| ) that Sun — ^ 0, 
which shows that S is closable. Similarly, T is closable if 5* is. O 

Let T and A be operators from X to Y, A being T bounded (|103|) . If 
T^^ exists and is a bounded operator from Y to X then AT~^ is an operator 
on Y to Y and is bounded by 

\\AT-^v\\ < a\\T-\\\ + b\\v\\ < (a||T-^|| + b)\\v\\. 

If 

a||T-^||+6<l (105) 

one can also show stability of bounded invertibility under relatively bounded 
perturbation. Indeed if ( p.05| ) is valid then S = T + A is automatically closed 
by Lemma ^. One has 



and thus 



S = T + A= {l + AT-')T. 
WT-H 



s-'\\ < 



l-a\\T 



-11 



In the Hilbert space one can show that selfadjointness is also stabile under 
relatively bounded perturbations. If T is selfadjoint and A is symmetric and 
T bounded with T-bound smaller than 1, then T + A is also selfadjoint. In 
particular T + A is selfadjoint if A is bounded and symmetric with D(A) D 
D(T) [|109|| . Note that the assumption that the bound be smaller than 1 
cannot be dropped in general. If T is unbounded and selfadjoint, and A = 
— T, then T-l-A is a proper restriction of the operator and is not selfadjoint 

An important example of a relatively bounded perturbation provides the 
Schrodinger operator in the 3-dimensional euclidean space R^ for a system 
of s particles interacting with each other by the Coulomb forces. In this case 
the formal Schrodinger operator is the 3s-dimensional Laplacian —A and the 
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perturbation V{x) has the form 



na:) = E? + E^, (106) 



where e, and e,i, are constants and 



'^j ~ l^3i-2 + -^sj-i + ^3j^ 



Tjk = [{x3j-2 - X3k~2y + {X3j-1 - X^k-l)'^ + {xsj - X^k 



,211/2 



It can be proved that the minimal operator T constructed from the formal 
operator —A is essentially self-adjoint with the self-adjoint closure Ho. If V 
denotes the maximal multiplication operator V{x) then it can be shown that 
V is relatively bounded with respect to Ho as well as to T with relative bound 
equal to zero 



4.1.2 Setting up the problem 

Henceforth we restrict ourselves to the ensuing problem. Let Ho be a closed 
linear operator acting in a Banach space and let V be another closed linear 
operator acting in this space, which is relatively bounded with respect to Ho- 
In addition let the bound satisfy ( |105| ) with T replaced by Ho- We shall be 
interested in the resolvent operator R\{z) := {Ho + XV — z)~^ defined for z 
not in the spectrum of Hx := Ho + \V (A G C). Usual perturbation theory 
starts from the identity 

Rx{z) = Ro{z)[l + XVRoiz)]-' (107) 

by developing the geometric series 

oo 

[1 + XVRoiz)]-' = Y.{-Xr[VRo{z)r. (108) 

n=o 

Under our assumptions, VRo{z) is a bounded operator, so this series will 
converge in norm for ||AV^-Ro(-2)|| < 1- If V^ is bounded then the condition 
can be replaced by ||Ay|| < d{z), where d{z) is the distance of z from the 
spectrum of Ho- 

Let us suppose that there exists a finite system S'(//o) of eigenvalues 
of Ho, separated from the rest of the spectrum of Ho by a closed curve 
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Figure 8: Perturbation of eigenvalues. 

encircling T!{Ho). Then there is a convex region A of the complex plane, 
containing the origin, such that for all A G A the spectrum T,{H\) of H\ = 
Ho + W is likewise separated by 6 into a part T,'{Hx) and a remainder. 
The eigenvalues in Tj'{H\) are analytic in A with only algebraic singularities 
[ |109|| . The usual Rayleigh-Schrodinger perturbation theory is obtained by 
substituting (|107[ - p.08| ) in the expression for Px, 

Px = ^J^Rxiz)dz, (109) 

the projection operator onto the subspace associated with the eigenvalues ex 
branching off from Eq in T,'{Ho), where F is a contour encircling ex but no 
other points of T,'{Hx) (see Fig. |). In the physical language if the contour 
F encloses all energy levels below the Fermi surface then Px is nothing but 
the Fermi projector or density matrix at zero temperature, respectively. 
To extend the region of convergence some summability method can be 



used | 108 , m j. Reeken used the Borel summability with f3 = 1. As we have 
shown, /i-method has a bigger region of convergence, so that an application 
of the method will provide a futher extension of the Rayleigh-Schrodinger 
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perturbation theory. Let A be a bounded operator acting in a Banach space. 
Consider an operator-valued function f{z), 



(1 + zA) 
Because of the Theorem ^ (see also 0), it is true that 

1 f^ , .^ i-zAt) 



/(^) = TT— ^ = / exp(-expt) ^ ^--^ dt (110) 

[1 + zA) Jo frl uin) 



for z G MLS{f). The r.h.s. of ( |110|) can be defined by means of the uniform 



convergence of Theorem |] as a hmit of entire functions. Due to Theorem 
for any compact set K C MLS{f) and each e there is to > and an integer 
N such that 



/(^)-E^tK/ e-^t-dt) 
±ri ain) Jo 



n=o /^W 

uniformly in z G K. Thus 

N 
— df = lim lim 

n=o K^) 



< e 



r e-'" y ^^^^ dt = lim lim Y i- z AY t° e-'" T dt 



as a limit of entire functions. Any operator valued function can be defined on 
its MLS in this manner. Thus this provides an alternative to the Dunford- 
Schwartz integral. To define an operator valued holomorphic function f{z) on 
MLS{f) as a limit of entire functions one can also use some other AMCSM 

The singularities of f{z) are just such z that —1/z is from the spectrum 
S(y4) of A. Therefore the r.h.s. of ( |110| ) will converge for z = 1 if there 
are no singularities on the segment [0,1], or, if —l/z G [— oo,— 1] ^ ^(^4), 
respectively. So, the following Corollary holds. 



Corollary 4 .■ Let A be a bounded operator acting in a Banach space. Let 
S(y4) be the spectrum of A. If S(y4) is contained in the complement of 
(— oo, —1], then 



AtY 



/■°° / XV- {-At " , 



'1 + ^) ■Jo ^o Kn) 
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Note that the domain of apphcation of the representation ( |110| ) strongly 
depends on the sign of the parameter A. 

Now, to extend the region of convergence of the Rayleigh-Schrodinger per- 
turbation theory one has to locate the spectrum of \VRo{z). The spectrum 
of XVRo{z) for z ^ E{Ho) is the set of all ^ = such that z G S[i7(-A//i)]; 



^ = may or may not belong to the spectrum ||111|| . To show this one uses 
identity 



[XVRoiz) - /i]-i = {-l/fi){H, - z)[H, + {-X/fi)V 



.1-1 



for /i 7^ 0. If 2; G S(/7(— A//i)) the second factor on the r.h.s. is not defined. 
Thus /i G T,{XV Ro{z)) . If z ^ Il(iJ(— A//i)) then the second factor is defined 
and its product with the unbounded operator Ho — z is a bounded operator 
on the whole space. 

From the above discussion and Corollary |^ we have the Lemma ^ which 
converts the problem of covergence of the /i-sum for a given A to the problem 
whether z belongs to the spectrum of some class of operators H^ or not. 



Lemma 9 .■ Let H^ and V be as above. Then 

[1 + XVRo{z)]~^ = / exp(-expt) ^[-A\/i?o(z)t]"//i(n)dt 

"'o n=o 

for all z which are not in the spectrum of H^ for v G {Xu | m G [0, 1]}. O. 

Note that if A is a real number, as usually happens in physical applica- 
tions, then one has only to consider a real parametric family of Hamiltonians 
in contrast to [|1 1 If . As a direct consequence of Lemma ^ and ( |109D one 



arrives at the following statement 
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Corollary 5 ; If for X G A no one of the branches emanating from Eq, since 
the perturbation parameter u varies on the straight line segment [0, A], crosses 
any other branches starting from other eigenvalues in Tj'{Ho), then a contour 
r exists encircling e^ and Sx but no other eigenvalues (see Fig. ^ such that 



P\= j exp(- exp t) ^ Rn 

Jo „_„ 



;-^*'"* 



/i(n) 
where Rn is the residue of Ro{z)VRo{z)...VRo{z) at Eo- <0 

4.2 Derivative analyticity relations 

Finally, we shall consider the problem of derivative analyticity relations 
(DAR) [|112| , |113|| . Let us briefly sketch the problem. The real and imag- 



inary parts of the forward scattering amplitude F{E) in high-energy physics 
(or in optics) are related by a dispersion relation of the form 

. ^ , 2E f°° Im F(E) , , , 

Re F(E)/E = - v.p. /^. ^^^^±4^iE, (111) 

where the pole terms and subtraction constants are for simplicity omitted. 
A major "shortcoming" of the dispersion relation is that one has to know 
the imaginary part Im F{E) on the whole infinite integration interval to 
obtain the real part Re F{E) at a given point. The problem of DAR is that 
of relating the real and imaginary part of F{E) at the same point. To cope 
with it in the general case is quite difficult. After the substitution x = lnE, 
Xo = lu-Eo, k = InE* and f{x) = Im F{E)/E the integral on the r.h.s. of 
(|111D takes the form, 



1 f°^ f{x + Xo) , 

- v.p. / ■ , . ^ dx. 

IT Jk-xo smh(a;j 

The solution of the problem of DAR is only known in a special case where 
f{x) can be analytically extended to an entire function ||113|, |114|| . Using 



the results of Section ^ one can establish the DAR for more general class 
of function /(x). In fact, under the assumption of regularity of f{z) on the 
real positive axis one arrives at the following representation of the dispersion 
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integral, 

1 /-OO roc ^ 

— v.p. / dx dt exp(— exp t) y^ / 

'1 «-' rC — Xq «-' O ,1-1 — r^ 



(n)| 



X. 



{xty 



n\fi{n) sinh(x) 



(112) 



because of Corollary |l[ If f{x) can be extended to an entire function, then 

J2'^=o^nf^^Kxo) converges, where 



l-OO ("OO 

/ dx (itexp(— expt) 

Jk—Xo Jo 



{xty 



n\fi(n) sinh(a;) 



One may then integrate term by term in ([L12|), reproducing the previous 



results ||113| , |114|| . For the references and recent status of the problem see also 
115|| . A different question is that of the practical use of the extension because 



the real positive axis is often the only place where a scattering amplitude is 
singular. 
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5 Conclusion 

We have found a family of moment constant summability methods fik which 
provide for k > 1 an analytic continuation of a function regular at the origin 
onto its whole Mittag-LefHer's star, in contrast to the Borel method (see 
Remark 2). These methods are intimately connected with the Borel one, 
since for a = P = 1 the Borel method is nothing but the /io-method. 

The methods discussed above can be succesfuUy applied to convergent 
as well as divergent perturbation series which diverge like (Inn)" or slower. 
One may encounter such a situation when regularizing a theory on a lattice 
[0, |117[ or using a cut-off in space time or momentum space 0, 0, |ll(j| ]. 



Of some interest is also an expansion in the "artificial" parameter proposed 
recently by Bender et. al. in the A04 theory and some other models [^] , 
and which seems to be convergent. When some additional information is 
at disposal or provided one deals with truncated series our method can be 
made more powerful by combineing it with a conformal mapping or the Fade 



approximation ||76|, |73, ^ |118|| . An interesting question is to compare a 



numerical efficiency of the /x^-methods with the Fade or some other methods 



[|119|| . So far we did not consider the possible generalization of the /ife-methods 



to more variables ||12(J|]. Of some interest is also the question whether the 



Wynn ^-algorithm for calculating the Borel integral could be adapted to 
the /ifc-methods [|121|| , and whether the inverse of a /i^-summable function is 



/ifc-summable [122 



It is obvious that the use of the above methods is not confined to per- 
turbation theory where the expansion parameter is a coupling. The method 
can be applied to the 1/A^ expansion ||5^ or to the e-expansion |^ as well. 



We have shown that the above given summability method solves the prob- 
lems (A) and (B), i.e., in the regular case provides an analytic continuation 
on the whole Mittag-Leffler (principal) star and can deal with the horn- 
shaped singularity as well. To our knowledge it is the first summability me- 
thod having these properties. We have discussed advantages and shortcomings 
of the method, applications to the Rayleigh-Schrodinger perturbation theory 
which cover the very important case of the Schrodinger operator in R^ for 
a system of s particles interacting with each other by the Coulomb forces 



108| , |109|| , and the derivative analyticity relations ||115|| . 

Like any analytic regular summability method our method may also have 
a wide domain of applications. We shall not give a list of them because the 
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reader can easily judge whether it is interesting for him or not. 

An open question still remains how the summability properties (both Bo- 
rel's and ours) are transported in equations such as the Dyson- Schwinger 
equations for Green's functions [|123|| and the problem of summing a pertur- 



bation series which violates SAC. To solve the last problem one has however 
to have more (nonperturbative) informations to pick up a physically plausible 
solution. 

One could also establish analogous results by other slightly modified 
Stieltjes moments of the form 



/i(n) := / exp(t — expt)t"(it , 

Jo 



etc. In this case one needs only to replace u^z) by uj{z) = exp z+\n{z+l) and 
Hr by Hji := {ReoJ{l/z) > ij{l/R)}. Note that in contrast to the previous 
case the boundary dH~^^ of H]^^ approaches the straight line Imz = 7i/2 
from below (see Fig. ^. 

We hope that a sufficient number of examples have been given to illustrate 
that to draw physical conclusions from the Borel (non) summability without 
knowing the region of analyticity to which a power series expansion should be 
summed up may sometimes be very dangerous. We have also shown that to 
prove the Callan-Symanzik assumption about the mass insertion terms one 
needs generalized SAC, derived in the paper, which apphes to horn-shaped 
regions. We then were able to formulate such conditions for a whole variety 
of horn-shaped regions. 

Finally, our method gives a generalization of [33] since Theorem |^ to- 
gether with Lemma ^ provide a summability mechanism which apart from 
invariance conditions and linear covariances preserves also nonlinear pertur- 
bative conditions such as the unitarity of the Feynman series. We note that 
Theorem ^ gives a generalization of SAC in regards to [^ ^ ^ . Indeed, 
by the theorem there cannot exist a function which 

i) is analytic in the horn H^ and continuous up to the boundary; 
ii) possesses there the asymptotic expansion ( |66D which has equal sign coef- 
ficients a„ for n > Ho- 

We would like to point out that assertions of the type that a quantity 
Qi equals another quantity Q2 to all orders of perturbation theory are very 
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vague unless the SAC are shown to be vahd. These SAC provide a complete 
generalization of Simon's work |Q. We have given arguments that four 
dimensional field theories (without UV cutoff) violate SAC and that probably 
there is no reason to look at the Borel transform of these thories. 



6 Acknowledgements 

First of all I should like to thank my parents for kind attention and sup- 
port. I should also like to acknowledged friutful discussions with members of 
our department. I am indebted to J. Fischer for continuous interest in my 
work and valuable comments, P. Hof ava for discussions on variety of physical 
problems, J. Chyla, P. Kolaf, S. Nemecek, and J. Rames for their help with 
computer facilities, and Mrs. M. Bouskova for drawing pictures. 

I take oportunity to thank J. Fuka for discussions about entire functions 
which help me in derivation of the properties of the new summability method 
in regular case as well as V. Sverak for many stimulating discussions which 
help me to derive the analogue of the Nevanlinna theorem for horn-shaped 
regions. Discussions with C. Klimcik on various aspects of field theories 
are also greatfully acknowledged. I thank J.Magnen, G.'t Hooft, and A.S. 
Wightman for encouraging correspondence, CPT, Ecole Polytechnique at 
Palaisseau for their warm hospitality and friutful discussions, E. Seller for 
kind hospitality in Miinich and discussions, as well as A. Smith Albion for 
reading a relevant part of the thesis. 



73 



7 Appendix 



Theorems 3.5.1-4 of the Phragmen-Lindelof type ||105|| can be easily modified 
to regions of the asymptotic type (1,77) (and also to other regions), i.e., if 
Fiz) 

i) is analytic in a region D of the asymptotic type (1, rj), and continuous 
in its closure D; 

ii)| F{z) \< Mexp[-S \ exp(l/r/2) |] z e D, 

S being an arbitrary positive constant, then F{z) = 0. 

So, if there exist two functions fi{z) and f2{z) satisfying the hypotheses 
of the theorem then their difference g{z) = fi{z) — f2{z) obeys the bound 

I g{z) |< ^p" I z I" exp[— exp w{n) + nlnw{n)] 

uniformly in z E D and n. Now, choosing n = [w^^ll/pr)] = [s{l/pr)] to 
optimize the bound, where [b] is now entire part of b and r =| z |, one gets 
that 

I g{z) |< Aexp[— exp(l/pr)] z E D, 

for region D of the asymptotic type {1,1]) with p < 1] . By virtue of ii) then 
g{z) = 0. 



On the other hand, whenever the bound ( pSf ) is relaxed, i.e., if in a domain 
D of the asymptotic type (1, r/) a bound 

I Rn{z) \<Ap^ I z \^ exp[- expw{N) + N\nw{N)] y z e D, 

is allowed, with p > rj, then for any p such that p > p > r] the function 
exp(— exp(l/p2)) satisfies the hypotheses of the theorem. In fact, 

sup I z~^exp(— exp(l/pz)) |=p^exp{ max Re\N In z — exp zW ~ 

z&D zel/pD 

~p^exp[-expw;(A^) + A^lnw(A^)]. 
However, this nontrivial function has trivial asymptotic series in D. <> 
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